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ABSTRACT 

The generation of a large-scale magnetic field in the kinematic regime in the absence of an a-effect 
is investigated by following two different approaches: the test-field method and the multiscale stability 
theory relying on the homogenisation technique. Our computations of the magnetic eddy diffusivity 
tensor of the parity-invariant flow IV of G.O. Roberts and the modified Taylor-Green flow confirm the 
findings of previous studies, and also explain some of their apparent contradictions. The two flows have 
large symmetry groups; this is used to considerably simplify the eddy diffusivity tensor. Finally, a new 
analytic result is presented: upon expressing the eddy diffusivity tensor in terms of solutions to auxiliary 
problems for the adjoint operator, we derive relations between the magnetic eddy diffusivity tensors that 
arise for mutually reverse small-scale flows v(x) and — v(x). 

Subject headings: MHD - magnetic fields - turbulence - dynamo 


1. Introduction 

It is well-known that at sufficiently high Reynolds 
number turbulence is characterised by a hierarchy of 
fluctuations interacting on a wide range of space and 
time scales. When this happens in a flow of conduct¬ 
ing fluid, magnetic field generation commences if the 
magnetic Reynolds number is sufficiently high (Mof- 
fatt 1978). As predicted by the magnetic induction 
equation governing the process of generation, small 
scales also develop in the generated magnetic field. 
The interaction of fine structures of flow and magnetic 
field usually influences the evolution of their large-scale 
parts. In particular, by Parker’s hypothesis, such an 
interaction may give rise to a mean electromotive force 
(e.m.f.), parallel to the large-scale magnetic field. 

In astrophysics, when the generation of the geomag¬ 
netic or solar magnetic field is under investigation, fine 
structures are generally of lesser interest than global 
ones. With present-day computers, it is impossible 
to resolve structures over the whole range of interact¬ 
ing scales; by choosing the domain of integration of 
the equations of magnetohydrodynamics, we can only 
focus on the large or small scales. However, in simu¬ 
lations of the global picture it is desirable to take into 
account the integral influence of physical processes at 
small scales. 


Since the 1960s, German scientists (Steenbeck, 
Krause & Radler 1966; see also Krause & Radler 1980) 
were developing the theory of mean-field electrodynam¬ 
ics (MFE), a first attempt supposed to advise how to 
do this. Perhaps, the best introduction to the ideas on 
which this theory is built is by one of its founders Karl- 
Heinz Radler (2007). The three-dimensional magnetic 
and flow velocity fields, b and v, are decomposed into 
“mean”, b and v, and “fluctuating”, b' and v', fields: 

b = b-|-b', v = v-|-v'. 

Any averaging procedure is deemed acceptable pro¬ 
vided it satisfies the Reynolds rules (see Radler 2007), 
e.g., planar averaging over any pair of Cartesian vari¬ 
ables, one-dimensional averaging along any given di¬ 
rection, or ensemble averaging for turbulent flows. The 
equations for mean magnetic field and fluctuations 
take the form 

_ _ _ 

= ryV^b -|-Vx(vxb-|-v'x b'), (1) 

at 

(9h' — 

= ryV^b' -f V x (v x b' -|- v' x b -|- (v' x b')'). 

^ ( 2 ) 

Here f' = f — f denotes the fluctuating part of a vector 
field f. The problem then reduces to the use of (2) for 
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expressing the mean e.m.f. v' x b' in terms of b and 
V. For simplicity, we henceforth assume that v = 0 
and v' is steady. In MFE, for homogeneous station¬ 
ary turbulence, the mean e.m.f. is usually expressed in 
terms of the mean magnetic held as 

v' X b'= yy (/C„(x-^,t-T)b(|,T) (3) 

-ICrf{x - I, t - r)V X b(|, r)) d| dr 

when averaging is planar (/Cq, and /C,, do not depend on 
the spatial variables over which the e.m.f. is averaged 
in the l.h.s.) — in general, rj should be dehned as a 
rank 3 tensor acting on Vb. Our task is to determine 
the kernels. In Fourier space, (3) implies 

-7^k.c;(v'X b') = a(k, u;)J'k,cdb - r 7 (k,a;)J'k,c;(V x b). 

(4) 

Here, following Brandenburg et al. (2008b), we have 
denoted 

= 11 f (x, t) dx dt, (5) 

Q!(k, w) = r7(k,a;) = 

In the limit k —> 0 and w —>■ 0, a and t] describe the 
(magnetic^) a-effect and eddy diffusivity correction^ 
tensors. 

The test-field method^ (TFM) for computing a and 
rj was developed within the MFE paradigm. To the 
best of our knowledge, it was hrst proposed by Schrin- 
ner et al. (2005, 2007). Perhaps, the most detailed 
description of the TFM procedure applied by Devlen 
et al. (2013) is found in Brandenburg et al. (2008a). 
The recipe is to solve equation (2) for zero-mean mag¬ 
netic perturbation b', where b is a test held. The ini¬ 
tial condition for b' can be any solenoidal small-scale 
zero-mean held (for instance, 0). For space-periodic 
magnetic helds, the test helds 

b = cos(k • x) e„ and b = sin(k • x) e„, (6) 

are chosen. By using sufficiently many independent 
test helds, we obtain a linear system of equations that 

^This paper is devoted to the study of magnetic o-effect and 
magnetic eddy diffusivity exclusively — as opposed to the hy¬ 
drodynamic a-effect known as the AKA-effect (see Frisch et al. 
1987; Dubrulle & Frisch 1991), or combined a-effect and eddy 
diffusivity emerging in large-scale perturbations of magnetohy¬ 
drodynamic regimes (see Chaps. 6-9 in Zheligovsky 2011). Note 
that the expression “magnetic a-effect” is sometimes used with 
a different meaning, designating a term proportional to current 
helicity that quenches against the kinetic a-effect. With this 
disclaimer in mind, we omit the attribute “magnetic” from now 
on when referring to the a-effect and eddy diffusivity. 

^We use here the terminology of the multiscale stability theory. 
In fact, the “corrections” can be much larger than the molecular 
diffusivity which they “correct” — the turbulent diffusivity can 
be by orders of magnitude larger than the molecular diffusivity. 

®Not to be confused with Kraichnan’s “test-field model” of turbu¬ 
lence (Kraichnan 1971), used by Sulem et al. (1975) as a method 
for closure of the hierarchy of moment equations. 


relates v' x b' through the unknown coefficients of a 
and T] to b. This system can be solved to obtain a 
and r). Similarly, the temporal dependence of the ker¬ 
nels in (3) can be “probed” in Fourier space by con¬ 
sidering the test fields 

b = cos(k • x) e”'^* e„ and b = sin(k • x) e”'^* e„. 

(7) 

In kinematic dynamo problems, where the evolu¬ 
tion of a weak magnetic held is studied (so that its 
inhuence on the flow via the Lorentz force can be ne¬ 
glected), the how velocity, v, is known a priori. It 
can be a stationary held, often supposed to have a 
vanishing average (v = 0), as have the hows that 
we consider in this paper. Alternatively, it can be a 
time-dependent how, for instance, supplied by an in¬ 
dependent hydrodynamic simulation. The kinematic 
dynamo problem is an instance of the full magneto¬ 
hydrodynamic (MHD) stability problem that focuses 
on the stability of non-magnetic states; the how and 
magnetic held perturbations then decouple since the 
Lorentz force is quadratic in the magnetic held. In a 
general setup, one considers the stability of an MHD 
regime featuring a non-vanishing magnetic held that 
affects the how, and therefore perturbations involve 
both the how and magnetic held that cannot be dis¬ 
entangled. 

MHD perturbations involving much larger spatial 
and temporal scales than those of the perturbed MHD 
regimes (which, e.g., can be periodic or quasi-periodic 
in space, and steady or periodic in time) can also be 
explored by an approach known as the multiscale sta¬ 
bility theory (MST). It originates from the studies of 
hydrodynamic stability (Dubrulle & Frisch 1991) and 
kinematic dynamo (Lanotte et al. 1999) and relies on 
mathematically precise asymptotic methods for ho¬ 
mogenisation of elliptic operators. An introduction to 
MST can be found in Zheligovsky (2011); the linear 
MHD stability problem for large-scale perturbations 
was considered by Zheligovsky (2003) (see also Chap. 6 
of Zheligovsky 2011). Here we will only consider the 
kinematic dynamo problem, and focus on the gener¬ 
ation of a magnetic held involving large scales by a 
small-scale huid how. For a steady how, the dynamo 
problem can be reduced to the eigenvalue problem for 
the magnetic induction operator: 

ryV^b -k V X (v X b) = Ab (8) 

(here rj denotes the magnetic molecular diffusivity and 
A is the eigenvalue). 

We assume that a large-scale magnetic mode 
b(X, x) depends on fast, x, and slow, X = ex, spatial 
variables, the how depends only on x, and the scale 
ratio s is small. We proceed by expanding a mode 
b(X, x) and the associated eigenvalue A (its real part 
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is the growth rate of the mode) in power series in £, 

OO OO 

b = ^b„(X,x)£”, A = ;^A„e”, (9) 

n—0 n—0 

and deriving a hierarchy of equations that the eigen¬ 
value equation yields in successive orders As it 
turns out, we can find each term of the expansions by 
solving successively equations from this hierarchy. For 
parity-invariant flows, that we will mostly consider, 
the series for the eigenvalue involves only even powers 
of e (see Section 3.5 of Zheligovsky 2011). 

The first equation in the hierarchy shows that the 
leading terms bp and Aq in the expansion (9) are, re¬ 
spectively, a small-scale eigenfunction and the associ¬ 
ated eigenvalue of the operator of magnetic induction. 
The asymptotic expansion can be developed for any 
eigenvalue Aq. For small scale ratios e, the growth 
rate may exceed Re(Ao) due to the interaction of the 
fluctuating components of the magnetic field and of 
the small-scale flow, but the corrections are at best 
linear in the small parameter s and hence small. We 
are mostly interested in the case where no small-scale 
magnetic field is generated and Aq = 0, since then 
the presence of large spatial scales can, in principle, 
result in the onset of magnetic field generation, i.e., 
in a qualitative change in the behaviour of the MHD 
system. (The case of an oscillatory small-scale kine¬ 
matic dynamo occurring for imaginary Aq was consid¬ 
ered in Section 3.8 of Zheligovsky 2011; it is, actu¬ 
ally, algebraically much simpler.) For Aq = 0, the first 
term bp is a linear combination of neutral small-scale 
magnetic modes with coefficients depending on the 
slow variable. These coefficients, called amplitudes, 
are determined from the solvability conditions for the 
higher-order small-scale equations from the hierarchy. 
When the problem is considered in a three-dimensional 
periodic domain, the kernel of the magnetic induc¬ 
tion operator comprises three neutral magnetic modes 
whose averages are the unit Cartesian coordinate vec¬ 
tors (generically, the kernel is three-dimensional). The 
amplitudes of these modes can clearly be interpreted as 
the Cartesian components of the mean magnetic field. 
Furthermore, by the theorem on the Fredholm alterna¬ 
tive (see, e.g.. Stone & Goldbart 2009), the solvability 
condition consists of the orthogonality of the inhomo¬ 
geneous term to the kernel of the operator adjoint to 
the operator of magnetic induction. Generically, this 
amounts to vanishing of the integral of the inhomo¬ 
geneous term over the periodicity box. As a result, 
when Aq = 0, equations for the amplitudes can be 
interpreted as mean-field equations, where the respec¬ 
tive terms describe the a-effect or the eddy diffusivity 
effect. 

The MST analysis reveals the non-universal charac¬ 
ter of (3). This asymptotic equality can be rigorously 
derived for a multiscale kinematic dynamo and vol¬ 


ume averaging in the generic case, when the kernel of 
the magnetic induction operator comprises three mag¬ 
netic modes with non-vanishing linearly independent 
averages. However, (3) does not necessarily hold for 
other types of averaging, or when the dimension of the 
kernel is higher — in the latter case, amplitudes of 
all neutral modes are involved in (3), as this happens, 
e.g., for translation-invariant convective dynamos (see, 
e.g., Chertovskih & Zheligovsky 2015). For MHD tur¬ 
bulence, (3) is likely to stem, for various averaging 
procedures, from the ergodic properties of the respec¬ 
tive MHD dynamical system, but, to the best of our 
knowledge, this equality was never fully demonstrated 
in the context of MFE at the mathematical level of 
rigour; it remains a phenomenological property of tur¬ 
bulence (such as, for instance, the Kolmogorov law). 

The standard a-effect and eddy diffusivity, arising 
in the limits k —>■ 0 and w —)■ 0, are an idealisation in 
which nonlinear terms, higher spatial derivatives and 
temporal derivatives of the magnetic field are omitted 
in the expression (3) for the e.m.f. This is justified if 
the mean fields vary sufficiently slowly in space and 
time, i.e., on scales much larger and longer than those 
of the fluctuations. While this simplification may be 
permissible in some cases, e.g., for forced turbulence 
with sufficient scale separation, for certain flows, such 
as the Roberts and Otani flows, it is not (Hubbard & 
Brandenburg 2009). A particularly striking example 
are flows H and HI of G.O. Roberts (1972); for de¬ 
scribing the nature of the dynamo in those flows, it is 
crucial to retain the convolution in time in the inte¬ 
gral operators in (3) (Rheinhardt et al. 2014). Then 
the electromotive force at a given time depends on the 
magnetic field also at earlier times, so the system pos¬ 
sesses “memory”. It is important to realise that the 
memory effect does occur even for steady flows such 
as those considered here. Excluding the memory ef¬ 
fect from consideration more often results in quantita¬ 
tive distortions, such as too high an estimate for the 
critical dynamo number (Rheinhardt & Brandenburg 
2012), rather than in qualitative changes. 

We note in passing that, instead of implementing 
an integral transform in both space and time, which 
is cumbersome, it is convenient to solve an evolution 
equation for the e.m.f. v' x b'. Such an equation 
was first derived by Blackman & Field (2002) using 
the T approximation, which captures temporal nonlo¬ 
cality, i.e., the memory effect (Hubbard & Branden¬ 
burg 2009). This was then extended by Rheinhardt 
& Brandenburg (2012) to capture also spatial nonlo¬ 
cality. Usually this also yields a satisfactory (at least 
qualitatively) description of the unusual phenomena 
related to the memory effect, such as the ones encoun¬ 
tered in flows 11 and HI of G.O. Roberts (Rheinhardt 
et al. 2014). These ideas will turn out to be important 
in Section 5.4, when we compare the magnetic field for 
the modified Taylor-Green flow (mTG) obtained from 
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direct numerical simulations (DNS) with that found 
by TFM. 

We have thus two independent theories: MFE, 
physical in spirit, especially when making simplify¬ 
ing assumptions regarding the kernel in the integral 
equation (3) for MHD turbulence, and MST, which 
yields a mathematically rigorous derivation of equa¬ 
tions for similar quantities from first principles. MST 
has a narrower scope, being applicable to treat only 
linear and weakly nonlinear MHD stability problems. 
While MST applies specifically to the limit k —> 0 and 
w —)■ 0, TFM can be applied to non-infinitesimal |k| 
and oj. It can therefore be used to assemble the kernels 
JCa and JCn- In other words, MST strives to describe 
an influence of the flow, characterised by certain tem¬ 
poral and spatial scales, on magnetic fields involving 
much larger scales; TFM is more ambitious in trying 
to assess the influence of both larger and smaller hy¬ 
drodynamic scales on magnetic field of a given scale. 
Although the limits k —)■ 0 and w —)■ 0 can be numer¬ 
ically expensive for TFM, a comparison with MST is 
possible. 

Recently Devlen et al. (2013) applied TFM to com¬ 
pute the magnetic eddy diffusivity in flows previously 
employed in the studies of Lanotte et al. (1999) and 
G.O. Roberts (1972) with the use of MST and a sim¬ 
ilar approach. Dubrulle et al. (2007) observed in sim¬ 
ulations the beginning of magnetic field generation by 
mTG when increasing the magnetic Reynolds num¬ 
ber starting from small values, which the authors cau¬ 
tiously attributed to the onset of the action of nega¬ 
tive magnetic eddy diffusivity investigated by Lanotte 
et al. (1999). Devlen et al. (2013) found, in agreement 
with G.O. Roberts (1972), that the so-called flow IV of 
G.O. Roberts (further referred to as R-IV) does yield 
negative magnetic eddy diffusivity, but they failed to 
reproduce the results of Lanotte et al. (1999) on the 
presence of negative magnetic eddy diffusivity in mTG. 
We resolve this controversy in the present paper and 
show that in a suitable parameter range eddy diffusiv¬ 
ity is negative, however, the relevant TFM averaging 
is not over the horizontal plane (which is applicable 
for R-IV), but one along the vertical direction, or a 
planar one over any of the other two Gartesian coordi¬ 
nate planes such that the average still depends on one 
of the two horizontal directions. 

The need for the cross-examination stems from the 
fact that some applications of the MFE ideas can fail to 
conform with the mathematical structure of problems 
under consideration. For instance, the mean e.m.f. 
computed as “an average over the lower half-volume, 
the upper half-volume, or, better still, one half of the 
difference of these two” was used in the studies of Gat- 
taneo & Hughes (2006, 2008) of the a-effect in con¬ 
vective dynamo in a layer. How could these proce¬ 
dures possibly help to track the evolution of the mean 


magnetic field? Such an averaging does not obey the 
Reynolds rules, namely, because averaging and taking 
the spatial gradient do not commute, and turns the 
midplane into an artificial boundary. In each half¬ 
cell the mean field depends only on the horizontal 
variables. The opposite a-effect values in two adja¬ 
cent half-cells force us to assume opposite mean fields 
over and below the midplane, in order to avoid sin¬ 
gularities in the a-effect operator at the midplane. 
This inevitably implies the existence of a boundary 
layer at the midplane. However, nothing resembling 
a boundary-layer kind of behaviour of magnetic field 
in the numerical solutions was reported ibid., clearly 
showing that averaging over a half-cell is unnatural 
and incompatible with the physics of the problem, 
and is also inappropriate from the mean-field electro¬ 
dynamics perspective. The a-effect operator must be 
calculated by averaging over the entire periodicity cell; 
the observed “antisymmetry of a about the midplane” 
(Gattaneo & Hughes 2006) simply implies that in these 
dynamos the relevant a-effect is zero (i.e., the a-effect 
operator is not involved in the equations for the evo¬ 
lution of the mean field), and the essential eddy effect 
is eddy diffusivity. Furthermore, the convective dy¬ 
namos considered ibid, are translation-invariant, and 
hence some amplitudes, essential in the description of 
the large-scale modulation of the generated instabil¬ 
ity modes, cannot be interpreted as mean fields^ (see 
Ghaps. 8 and 9 in Zheligovsky 2011); neglecting these 
modes is also likely to affect the results of Gattaneo 
& Hughes (2006, 2008). As a result, no sound conclu¬ 
sions concerning the a-effect, intended for astrophysi- 
cal or general MHD applications, can be drawn from 
the findings of those two papers. 

Our paper is organised as follows. In Section 2 
we remind the reader of the MST formalism for the 
large-scale kinematic dynamo. In Section 3 we calcu¬ 
late, in the MST framework, the operator of magnetic 
eddy diffusivity for R-IV using its many symmetries, 
and state results of the computation of its two coef¬ 
ficients. In Section 4 we discuss how the symmetries 
of mTG reduce the number of auxiliary problems in¬ 
volved in MST computations of eddy diffusivity, and 
present numerical results. Despite using algorithms 
that differ drastically from those used by Lanotte et 
al. (1999), we reproduce the results of this paper with 


^We will see in Section 2 that a large-scale magnetic mode has the 
structure b = B"(X)S„(x)-l-Ofe), where N = dimkeril is 

the number of independent small-scale neutral magnetic modes 
S„(x). By normalizing the small-scale modes, we can impose 
the conditions 

e„ for n < K, 

0 for iy -I- 1 < n < iV. 

We then find b = B"(X)en; thus, for n < K the am¬ 
plitudes i?"(X) have the sense of the mean components of the 
mean field b; for n > K + 1 no any such or similar interpretation 
is possible. 
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4 significant digits. In Section 4.2 we explain why no 
large-scale dynamo was found for mTG by Devlen et 
al. (2013), and show that eddy diffusivities obtained 
by TFM with an alternative planar averaging quali¬ 
tatively agree with the MST values. In Section 4.3 
we show that the growth rates of large-scale dynamo 
modes have the symmetry properties implied by the 
structure of the eddy diffusivity operator. In Section 5 
we demonstrate that the TFM procedure with the spa¬ 
tial averaging reproduces the MST a-effect and eddy 
diffusivity tensors, and consider analytically and nu¬ 
merically the difference of the two approaches for a 
planar averaging using mTG as an example. Gonclud- 
ing remarks end the paper. 

2. The mathematical theory of generation 
of large-scale magnetic field 

We review here the results of application of MST 
for the investigation of large-scale magnetic field gen¬ 
eration by small-scale steady flow of electrically con¬ 
ducting incompressible fluid (Lanotte et al. 1999; Zhe- 
ligovsky et al. 2001; Zheligovsky 2011). We consider 
the kinematic dynamo problem as a problem of deter¬ 
mination of the spectrum of the magnetic induction 
operator, which enables us to hnd growing large-scale 
modes even when in addition a small-scale dynamo op¬ 
erates. For the sake of simplicity, both the large-scale 
magnetic mode b(X, x) and the flow v(x) are assumed 
to be 27r-periodic in each fast spatial variable Xi. The 
mode is solenoidal and satisfies the eigenvalue equa¬ 
tion (8) for the magnetic induction operator. 

1. Magnetic a-ejfect. Generically, the average of 
the leading term in the expansion (9) of a magnetic 
mode, B(X) = (bo(X,x)), and the leading term in 
the expansion of the associated eigenvalue, A = Ai, 
are a solution to the eigenvalue problem for the a- 
effect operator, 


is the small-scale magnetic induction operator, and e„ 
are unit vectors of the Gartesian coordinate system. 
S„(x) are solenoidal. 

Let B(X) be a solenoidal space-periodic solution to 
the eigenvalue problem (10) whose associated eigen¬ 
value is A. Then B(/iX) is also a solenoidal solution 
to (10) whose associated eigenvalue is pA; for any in¬ 
teger fj,, positive or negative, this mode possesses the 
spatial periodicity of the original mode B(X). Thus, 
a mean field, that is initially an infinite sum of modes 
defined by (10), grows in general superexponentially; 
consequently, the large-scale magnetic field grows and 
destabilises the MHD system on time scales that are 
intermediate between the fast time t and the slow time 
T = et (unless all modes defined by (10) are associated 
with imaginary eigenvalues A). 

2. Magnetic eddy diffusivity. A field f is parity- 
invariant, if 

f(-x) = -f(x), (13) 

and parity-antiinvariant, if 

f(-x)=f(x). 

For parity-invariant flows v, parity-invariant and 
parity-antiinvariant vector fields constitute invari¬ 
ant subspaces of the magnetic induction operator 
Hence, vector fields S„(x) are parity-antiinvariant, 
and the a-effect is absent: 21 = 0. The magnetic 
field (9) is then 

b(X,x) - ;^(H"(X)(S„(x)+e„) (14) 

n^l \ 

(X)G™„(x)^ + 0(£2), 

where vector fields Gmn(^) are zero-mean solutions to 
auxiliary problems of type II: 




m—l 




Vx X 21B = AB, (10) 

in the subspace of solenoidal fields, Vx • B = 0. Here 
the tensor of magnetic a-effect, 21, is the 3x3 matrix 
whose n-th column is (v x S„), (•) denotes the aver¬ 
age over the periodicity cell = [0, 27r]^ of the fast 
variables, 


(f)(X) = (27r)-3 / f(X,x)dx, 

vector fields S„(x) are zero-mean solutions to auxiliary 
problems of type I: 


-CG™„ = -277—At - e™ X (v X (S„ -h e„)). (15) 

OXm 

Gmn(x) are parity-invariant. 

The solenoidal mean part of the leading term in the 
expansion (9) of the mode, and the leading term in the 
expansion of the associated eigenvalue, A = A 2 , are a 
solution to the eigenvalue problem for the operator of 
magnetic eddy diffusivity: 


ryVxB -I- Vx x ^ ^ ) S)r 


n—1 m—l 


dB'^‘ 


= AB. 


(16) 




£S„ = - 


dw 


dxn 

ii(S„ -I- e„) = 0, 


( 11 ) 

( 12 ) 


Here, 3? is the tensor of eddy diffusivity correction, 
^mn = (v X Gm„). (17) 


iib = TyV^b -I- Vx X (v X b) 
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We assume that the mean fields reside and are 
bounded in the entire space Hence, solutions to 



the eigenvalue problem (16) are Fourier harmonics^ 

B(X) =Be‘'i'^, B-q = 0. (18) 

Here, B = and q = ( 91 ,( 72 , 93 ) are con¬ 

stant vectors satisfying the conditions |q| = 1, B-q = 0 
(solenoidality of the mean magnetic mode) and 

3 3 

-ryB - q X X! X! ®mnS”9m = AB. (19) 

n—1 m—1 

Solenoidality of the modes implies 

B = PtT + /3pP, ( 20 ) 

where 

T = (- 92 , 91 , 0 ), P = (9193, 9293 ,-(9?+92)) ( 21 ) 

(this is equivalent to decomposing the mode into the 
toroidal and poloidal (this is equivalent to decompos¬ 
ing the mode into toroidal and poloidal components). 
Substituting (20) into (19) and scalar multiplying by 
T and P, we recast (19) into an equivalent eigenvalue 
problem in the coefficients Pt and Ppi 

- Y, + PpP^)qm = ( 9 ? + 9l)(ry + A)Pu 

m,l,n 

( 22 ) 

^ DL„T'(AT" + /?pP") 9 ™ = (9? + ql){v + A)/ 3 p. 

m,l,n 

(23) 

Taking into account the symmetries of the generat¬ 
ing flow can considerably simplify the eigenvalue prob¬ 
lem (22)-(23) (see Sections 3 and 4). 

Eigenvalues A depend on the wave vector q of the 
large-scale amplitude modulation: A = A(q). If the 
real part of A(q) is the maximum of Re(A(q)) over unit 
wave vectors q, then rj^ddy = —A(q) is called the min¬ 
imum magnetic eddy diffusivity. When Re(r7eddy) > 0, 
generation of large-scale magnetic field by the mech¬ 
anism of negative eddy diffusivity is possible. From 
a physicist’s point of view, this mechanism is impor¬ 
tant only if the flow v does not generate small-scale 
magnetic fields (i.e., fields of the same spatial period¬ 
icity, as that of the flow), because otherwise small-scale 
magnetic fields grow and destabilise the MHD system 
on time scales of the order of unity, which is faster 
than the growth of the large-scale field in the slow 
time T = e^t. This can also be interpreted as follows: 
when only the small-scale dynamo is acting, the mag¬ 
netic field can involve Fourier harmonics of arbitrarily 
large wave lengths (compatible with the boundary con¬ 
ditions, i.e., not exceeding the size of the periodicity 


®The vector eq is analogous to the wave vector k referred to in 
the exposition of TFM in the Introduction. 


box when periodicity conditions in space are consid¬ 
ered) , but they decay and are unimportant for genera¬ 
tion. By contrast, when the small-scale dynamo is in¬ 
active, the presence of large scales in the field becomes 
a key ingredient, without which the mechanism of neg¬ 
ative eddy diffusivity cannot make a dynamo work. It 
can also happen that the small- and large-scale mech¬ 
anisms coexist and are acting simultaneously. 

3. Computation of the eddy diffusivity tensor. The 
load of computation of the tensor of eddy diffusivity 
correction is halved, if instead of computing the fields 
Gmn one solves auxiliary problems for the adjoint op¬ 
erator (Zheligovsky 2011): 

= V X ep (24) 

for zero-mean fields Z/, 1 < 7 < 3, the adjoint operator 
being 

£* z = rjV^'L — V X (Vx X z), 
since, as it is easy to see from (17), (15) and (24), 

= (^i ■ -b X (v X (S„-b e„))^^ . 

(25) 

4 . Relations between tensors of magnetic eddy diffu¬ 
sivity correction for mutually opposite flows. The av¬ 
erage (25) can be expressed in terms of solutions to the 
auxiliary problems for the adjoint operator. We dec¬ 
orate by the superscript “minus” the quantities perti¬ 
nent to the reverse flow —v: 

= ryV^b — Vx x (v x b), 

£-(S-+e„) = 0, (ii-)*(Zr+e0=0. 

Clearly, (24) implies 

(Vx X Z/ -b e;) = 0, (26) 

and hence for all I, 

VxxZz = Sf ^ Zi=p-^WZ^{vx{Sf +ei)), 

(27) 

where denotes the inverse Laplacian in the fast 
variables. Using the analogues of these relations for 
the flow V to eliminate S„ in (25), we obtain 

S)U = ^(z,.(2VxX ^-e™x V^Z-^y (28) 

Applying standard vector analysis transformations, we 
can express this average as an integral of the scalar 
product of Z~ and a field resulting from the action of 
a differential operator on Z;. By self-adjointness of the 
Laplacian and the curl, and antisymmetry of the triple 
product with respect to permutation of its factors, we 
find 

= -(®-)L- (29) 
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When small-scale magnetic fields are not generated 
(i.e., all eigenvalues of the small-scale magnetic induc¬ 
tion operator have non-positive real parts), the auxil¬ 
iary problems can be solved numerically by comput¬ 
ing S„ -|- e„ and Vx x Z; -|- e; as small-scale dominant 
eigenmodes of the magnetic induction operators XI and 
Xl~, respectively, (see (12) and (26)) in the subspace 
of solenoidal vector fields whose average can be non¬ 
zero. The same small-scale eigenvalue code is applied 
to solve all these six eigenproblems, the flow being re¬ 
versed, V — >■ —V, when computing Vx x Z;. 

3. Generation of large-scale magnetic field 
by R-IV 

G.O. Roberts (1972) studied how simple flows de¬ 
pending on two spatial variables xi and X 2 (deemed 
horizontal), such as (30) (see below), generate mag¬ 
netic fields, whose dependence on x^ enters via the 
factor Here, e is a small parameter; thus this 

work is clearly in the multiscale spirit, although he 
did not present the complete multiscale formalism, nor 
derived the operator of eddy diffusivity. His flow IV 
(labelled here R-IV) lacks the a-effect; it is the first 
known example of a dynamo exploiting the mechanism 
of negative eddy diffusivity, as was suggested previ¬ 
ously on general grounds (Zheligovsky et al. 2001). To 
the best of our knowledge, Devlen et al. (2013) were 
the first to identify and study in detail this mechanism 
for R-IV. It should be emphasised that flows H and HI 
are also non-helical dynamos, thus indicative of a neg¬ 
ative eddy diffusivity effect; however, later those flows 
turned out to have positive eddy diffusivity, and their 
dynamo action was identified as being due to turbulent 
pumping with a time delay (Rheinhardt et al. 2014). 

We follow Devlen et al. (2013) in investigating large- 
scale generation by R-IV. In the spatial variables in¬ 
troduced by Tilgner (2004) (rotated by 45° about the 
vertical axis with respect to the variables used by 
G.O. Roberts 1972), its Cartesian components are 

vi = Vasina:! cosx2, 

V2 = —v^cosa;! sina;2, (30) 

V 3 = sinxi. 

It is clearly incompressible and parity-invariant (see 
(13)), thus lacking an a-effect. 

3.1. The effect of symmetries 

The symmetries of the flow control the structure of 
the tensor of eddy diffusivity correction 3D. 

1. Translation antiinvariance with respect to the 
shift by half a period in xi of R-IV: 

v(a;i, a;2,3:3) = -v(a;i -b tt, X2, X3). 

(Note that the nonlinearity in the Navier-Stokes equa¬ 
tion is not invariant for the antisymmetry of this 


type, making this choice of flow somewhat academic.) 
Hence, applying the operation of shift by half a period 
in the direction xi, which we denote by V 

f(a;i,a;2,a;3) = f(a;i -b tt, 3:2,0:3), 

to the eigenvalue problem (26), we find 

Z- = Z„. (31) 

Substituting this into (28), using the self-adjointness of 
the Laplacian, the curl and operator and integrating 
by parts in Xm the first term in (28), we obtain 

(32) 

and = 0 for any flow possessing translation anti¬ 
invariance with respect to the shift by half a period in 
one of the spatial variables. 

2. Symmetry in X 2 of R-IV: 

V^{xi, -X2,Xz) = v'^{xi,X2,xz), 

v‘^{xi, -X2 ,xz) = -v‘^{xi,X2,x:i), (33) 

v^(a;i, -a; 2 ,3:3) = v^{xi,X 2 , X3); 

antisymmetry in 3 : 2 , is defined by changing here the 
signs in the r.h.s. to the opposite ones. Clearly, the 
curl or vector multiplication by R-IV maps fields, sym¬ 
metric in 3 : 2 , to fields, antisymmetric in X 2 , and vice 
versa. Consequently, fields symmetric and antisym¬ 
metric in X 2 constitute invariant subspaces of the op¬ 
erators of magnetic induction Xi and Xi“. It follows 
from (11) that S„ are symmetric in X 2 for odd n and 
antisymmetric in 3:2 for n = 2; (24) implies that Z; are 
antisymmetric in X 2 for odd I and symmetric in X 2 for 
I = 2 . 

Vector multiplication by also maps symmetric 
in X 2 fields to antisymmetric ones and vice versa for 
odd m, and does not change the symmetry and anti¬ 
symmetry of a field in X 2 for m = 2. Therefore, (25) 
implies 

®rnn = 0) if ? -b TO -b n is odd. (34) 

3. Wave vector parity. We call “even” a three- 
dimensional vector field depending on two spatial vari¬ 
ables xi and 32, when it is a linear combination of har¬ 
monics Bq e“^°° such that = 0 if gi -b 92 is even and 

= 0 if -b <72 is odd; we call a fieldJ‘odd”, 
when it is a linear combination of harmonics Bq e“^°° 
such that R^ = 0 if <71 -b <72 is odd and R^ = R^ = 0 if 
<71 + <72 is even. Clearly, in this terminology R-IV (30) 
is even. 

Taking the curl or calculating the vector product 
with R-IV transforms an even field into an odd one, 
and vice versa. Thus, even and odd fields constitute 
invariant subspaces of the magnetic induction opera¬ 
tors XI and Xl~. By virtue of (11) and (24), S„ are even 
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for n = 1, 2 and odd for n = 3, while Z; are odd for 
I = 1,2 and even for I = 3. Vector multiplication by 
maps odd fields into even ones and vice versa for 
m = 1 , 2 , and it does not change this type of “parity” 
for TO = 3. Using this, it is easy to show that 

®11 = = 0. (35) 

4- Swapping of the horizontal coordinates xi 0 x 2 . 
Since the flow and solutions Si and S 2 to auxiliary 
problems of type I are independent of the vertical coor¬ 
dinate, equations for horizontal components of Si and 
S 2 involve the vertical components neither of the flow, 
nor of the respective S„. We establish by inspection 
that the held (S'|(a: 2 ,a;i + tt), S\{x 2 ,xi -f tt)) satishes 
the same equation as (S'! (x). S'! (x)), and hence 

Sl{xi,X2) = Sl{x2,Xi -I- Tt), 

Sl{xi,X 2 ) = Sl{x 2 ,xi + tt). (36) 

We use the second of these relations to show that 

(37) 

Denote ip = v^sinxi sinx 2 ; clearly, 

(u\u^,0) = Vx X (f/'ea). 

Since the how is independent of X 3 , for n = 3 the 
source term in the r.h.s. of ( 11 ) vanishes, and hence 
S^ = 0. Therefore, by virtue of (27) for I = 3, 
Z 3 = ( 2 r 7 )-iV xV’- Since gradients are orthogonal 
to solenoidal helds in the Lebesgue space, in expres¬ 
sion (25) for 2)21 the term involving the derivative 
dSi/dx 2 is zero. Hence, on the one hand, 

2)21 =(277)“^ (Vx X (■i/'e2) • (v X (Si -f ei))) 
=(277)“^ {ipe2 • Vx X (v X (Si -f ei))) 

= - (2r7)"VV'e2 ■ 

We have used here the self-adjointness of the curl, (12) 
for n = 1 and the self-adjointness of the Laplacian. On 
the other, by the self-adjointness of the curl and by 
virtue of (25), (31) for I = 2, (12) for n = 2 and the 
relation S 3 = 0 , 

2)?3 = (Z2 • (ei X (v X 63))) 

= - (Vx^(v X (§2 -f 62)) • (ei X Vxt/’)^ 

=^“'(Vx"(^V2S2)-7/-ei) 

Thus, (37) follows from (36). 

5. Eddy diffusivity. We calculate now eigenvalues of 
the eddy diffusivity operator (16). By (32), the sums 


involving ftp and j3t in the l.h.s. of (22) and (23), re¬ 
spectively, vanish, and therefore these equations yield 
the same eigenvalue 

A = -77 + (g? + g|)-i ^ 2)(„„T'P"g„. 

m,Z,n 

By virtue of (21), (32), (34), (35) and (37), 

A = -77 -f D?2 (9? + dl ) + ® 3 i <&■ 

Actually, we have calculated the symbol of the eddy 
diffusivity operator acting on mean fields (defined by 
the l.h.s. of (16)); hence this operator for R-IV (30) is 

+ +(^-®3i)^- (38) 

The minimum eddy diffusivity is 

?7eddy = ?7 + niin(-2)?2, -©31)- 

3.2. Numerical results 

The coefficients 77 — 2)32 and — of the eddy dif¬ 
fusivity operator (38) have been computed using (25). 
Solutions S„ to auxiliary problems of type I, and solu¬ 
tions Z; to auxiliary problems for the adjoint operator 
have been computed by optimised iterations (Zheligov- 
sky 1993) as the dominant (associated with the zero 
eigenvalue) eigenfunctions of the operators of magnetic 
induction 21 (12) and 21~ (26). Iterations were termi¬ 
nated when the estimate of the dominant eigenvalue 
was below 10“^*^ in absolute value and the norm of 
the discrepancy for the normalised associated eigen¬ 
vector was below 5 • 10“^^. A resolution of 64^ Fourier 
harmonics was used before dealiasing, that was per¬ 
formed by discarding harmonics with wave numbers 
over 28. With this resolution, energy spectra of solu¬ 
tions to auxiliary problems decay by 30 orders of mag¬ 
nitude for 77 = 0.2 and still by 4 orders for 77 = 0.01. 
Plots of 77 — 2)32 and 77 — 2)33 are shown in Fig. 1 for 
0.01 < 77 < 0 . 2 ; in this range of molecular diffusivi- 
ties no generation of small-scale magnetic fields takes 
place. 

Figure 1 implies that a large-scale magnetic field 
is not generated for horizontal wave vectors q of the 
harmonic large-scale modulation, but it is generated 
for the vertical wave vector. We did not check if gen¬ 
eration of small-scale fields starts on further decreas¬ 
ing the molecular diffusivity; the behaviour of plots 
in Fig. 1 suggests that it may take place, and then 
—>■ — 00 and 77 — 2)33 —>■ 00 when 77 approaches 
the critical value for the onset of small-scale magnetic 
field generation. If this happens, the type of the gener¬ 
ated large-scale field changes: for smaller 77, generation 
of large-scale magnetic field for the vertical wave vec¬ 
tor q replaces the one for horizontal wave vectors. 

Plots of the two entries, 77 — 2)^2 and V — ©li, of the 
eddy diffusivity tensor are shown in Fig. 1 for a range 





Figure 1. Entries 77 — £>22 (upper panel) and rj — 
(lower panel) of the eddy diffusivity operator (38) for 
R-IV (30). 

of molecular diffusivities over the critical value for the 
onset of generation of the small-scale magnetic field. 
The form (38) of the operator of eddy diffusivity cor¬ 
roborates the conclusions of Devlen et al. (2013) that 
the eddy diffusivity tensor for R-IV is diagonal and has 
a double eigenvalue, i.e., its action on fields depend¬ 
ing on the vertical slow variable (which was the object 
of the studies of G.O. Roberts 1972 and Devlen et al. 
2013) is homogeneous. However, since the two coeffi¬ 
cients in (38) are distinct (see Fig. 1), eddy diffusivity 
is anisotropic, differing in the vertical and horizontal 
directions. Comparison of the lower panel of Fig. 3 
in Devlen et al. (2013) with the plot of ry — in 
the lower panel of Fig. 1 reveals a reasonable qualita¬ 
tive consistency between the large-scale magnetic field 
growth rates, obtained by Devlen et al. (2013) in DNS, 
and the MST minimum eddy diffusivity values, shown 
in the lower panel of Fig. 1, for roughly ry < 0.1. How¬ 
ever, while here we study eddy diffusivity in the limit 
e —> 0, Devlen et al. (2013) computed turbulent mag¬ 
netic diffusivity for finite scale separations; in partic¬ 
ular, the plot in Fig. 3 ibid, shows the diagonal entry 
of ri{e) for £ = 1 , whose behaviour is clearly different 


from that of ry — presented in Fig. 1. 

4. Generation of large-scale magnetic field 
by the modified Taylor Green flow 

As Lanotte et al. (1999) and Devlen et al. (2013), 
we now consider large-scale generation by the modified 
Taylor-Green flow (mTG), whose components are 

vi = sin xi cos X 2 cos X 3 + a sin 2 x 1 cos 2 x 3 

-I- b cos X 3 (sin xi cos 3 x 2 + c sin 3xi cos X 2 ), 

V 2 = — cos xi sin X 2 cos X 3 -I- a sin 2 x 2 cos 2 x 3 (39) 

— &cosx3(cos3xi sinx 2 -|- ccosxi sin 3 x 2 ), 

V 3 = — a sin 2 x 3 (cos 2 xi + cos 2 x 2 ) 

-I- d sin X 3 (cos Xi cos 3 x 2 ~ cos 3xi cos X 2 ). 

The flow is incompressible for d = b{3c— 1), which will 
be henceforth assumed. We now consider its symme¬ 
tries relevant for simplification of the eigenvalue prob¬ 
lem (22)-(23) and calculate the eigenvalues. 

4.1. The effect of symmetries 

1. Symmetries in Xi. A field f = (/^,/^,/^) is 
called symmetric in Xi, if for all i and j such that 

l<bj<3 

fii-lfx^, {-lfx2, (-!)'• X3) = (-!)'• /^'(x) 

(cf. (33)), and antisymmetric in Xi, if for all such i 
and j 

f{{-lf'Xi,{-lY^X2,{-lf'X3) = (-l)^"*^* /^(x), 

where is the Kronecker symbol. Since mTG is sym¬ 
metric in all Xi, it is parity-invariant and lacks an 
a-effect. 

When a flow is symmetric in Xi, vector fields pos¬ 
sessing the symmetry or antisymmetry in Xi constitute 
invariant subspaces of the operators of magnetic induc¬ 
tion £ and . Since all the three symmetries in Xi are 
independent, there are eight such invariant subspaces. 
We label them by 3-character strings; A and S in the 
7 -th entry of the label indicate that vector fields in the 
invariant subspace are symmetric or antisymmetric in 
Xi, respectively. For instance, SAA labels the invariant 
subspace, in which vector fields are symmetric in xi 
and antisymmetric in X 2 and X 3 . 

By virtue of (11) and (24), invariance of the fields, 
symmetric or antisymmetric in Xi implies that S„ for 
n i and Zi are symmetric in Xi, while Si and Z; for 
I ^ i are antisymmetric in Xi. Gonsequently, = 0 
when none of the indices I, n and m are equal to i. It 
follows 

=0 if m = n, OT I = m, OT I = n. (40) 
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Figure 2 . An elementary periodicity cell of mTG ( 39 ): 
a prism whose edges are periodicity vectors Ci ( 42 ). 
The vertex O' of the upper square base O'A'B'C' 
projects down along the vertical into the centre Q of 
the lower square base OABC of the prism. 

A Swapping of the horizontal coordinates xi ^ X2- 
The mTG has also a symmetry, which we denote by 7: 
a field f is 7-symmetric, if 

f^ixi,X2,X3) = f^ix2,Xi,X3 + Tt), 
f^ixi,X2,X3) = f^ix2,Xi,X3 + Tt), 
f^ixi,X2,X3) = f^ix2,Xi,X3 + Tt), 

and 7-antisymmetric, if 

f^{xi,X 2 ,X 3 ) = -/^(x 2 ,a;i,a ;3 + tt), 
f‘^{xi,X 2 ,X 3 ) = -f^ix 2 ,Xi,X 3 + 7 t) , 
f{xi,X2,X3) = -f{x2,Xi,X3 + 7 t) . 

7-symmetric and 7-antisymmetric fields constitute in¬ 
variant subspaces of the operators of magnetic induc¬ 
tion S, and ST. This implies that S3 is 7-symmetric, 
Z3 is 7-antisymmetric, for n = 1,2 the field S„ 
is mapped by the 7-symmetry to S3_„, and Z„ is 
mapped by the 7-antisymmetry to Z3_„. (We thus 
need to compute just 4 solutions to the auxiliary 
problems, say, Si, S3, Zi and Z3; S2 and Z2 can 
then be obtained by applying the 7-symmetry and 
7-antisymmetry to Si and Zi, respectively.) Conse¬ 
quently, the remaining non-zero entries of the eddy 
diffusivity correction tensor satisfy the relations 

®12 = ®23 “ ®31 “ ~® 32 - ( 41 ) 

When the 7-symmetry acts on a vector field, the 
symmetry or antisymmetry in xi becomes a symmetry 
or antisymmetry, respectively, in X2, and vice versa. 
Thus, the 7-symmetry maps ASA and SAA mutually one 


into another, as well as ASS and SAS. Since it also maps 
an eigenfunction of the operator of magnetic induction, 
XI, to an eigenfunction, restrictions of Xi on the two 
invariant subspaces, constituting any of the two pairs, 
have the same spectra. The subspaces AAA, AAS, SSA 
and SSS are invariant under the action of the symmetry 
7 ; each of them splits into invariant subspaces of Xi, 
that consist of 7 -symmetric or 7 -antisymmetric fields. 

3. Wave number parity. Inspection of (39) reveals, 
that mTG is comprised of Fourier harmonics e'*^'^ in 
which all the three wave numbers ki have the same 
parity, e.g., the sum of any two wave numbers is even. 
Consequently, the obvious periodicity cell of the 
flow, which is a cube of size 27r whose edges are parallel 
to the Cartesian coordinate axes, is not the smallest 
one. It is easily seen that a flow possessing the parity 
property of this kind is invariant under shifts along 
any of the periodicity vectors 

Cl = (7r,7r,0), C2 = (7r,-7r,0), Cs = (tt, 0, tt). 

.(42) 

(Clearly, this translation invariance implies 27r-periodi- 
city in any Cartesian variable Xi.) Therefore, elemen¬ 
tary periodicity cells of the flow are prisms whose edges 
are these vectors (see Fig. 2). Alternatively, one can 
regard the parallelepiped 

0 < Xl < 27r, 0 < X2 < TT, 0 < 0:3 < TT 

as an elementary periodicity cell of the flow, assuming 
the “brick wall” tiling of space by these cells, in which 
the parallelepipeds are arranged in infinite “bars” par¬ 
allel to the xi-axis, and any two adjacent bars are 
shifted along the a;i-axis by half a period relative each 
other. The volume of the elementary periodicity cells 
of both types is 27r^, e.g., a quarter of that of T^. Nev¬ 
ertheless, by a small-scale dynamo we understand the 
generation of magnetic fields which are 27r-periodic in 
each variable Xi. 

Each invariant subspace of Xi considered above fur¬ 
ther splits into subspaces of the so-called even and odd 
fields that are linear combinations of Fourier harmon¬ 
ics such that the sums of the wave numbers ki + kj are 
even or odd. We therefore extend the labels of invari¬ 
ant subspaces by two additional characters denoting 
the parity of the sums ki k2 of wave numbers in 
the horizontal directions (the fourth character), and 
the sums ki -\- ks of wave numbers in directions xi 
and X3 (the fifth character); E and 0 indicate even or 
odd such sums, respectively. For instance, the invari¬ 
ant subspace SAAOE consists of vector fields that are 
symmetric in Xi, antisymmetric in X2 and X3, and are 
comprised of Fourier harmonics such that the sum of 
wave numbers in the horizontal directions is odd and 
the sum fci -|- ^3 is even; the spectrum of £ is the same 
in this subspace and in ASAOO. 

4 . Eddy diffusivity. For an eddy diffusivity cor¬ 
rection tensor with the properties (40) and (41) stem- 
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ming from the symmetries of the flow, in Xi and 7 , it 
is straightforward, using (21), to reduce (22)-(23) to 

(-7 + + ql) + '^hqDPt = 

(-77 + £> 23(91 + 92 ) + '^IiqDPp = A/3p, 

respectively. Therefore, the two eigenvalues are 

Ai = —77 + £^ 2(91 + 92 ) + ®3i93> (43) 

A 2 =-77 + £ 23 ( 9 ? + ^ 2 ) + ®3i93 (44) 

(this explains Fig. 3 in Lanotte et al. 1999). The min¬ 
imum of — Ai(q) and — A 2 (q) over unit wave vectors 
q occurs for the vertical unit vector q = 03 or at any 
horizontal unit vector q = (< 71 , ( 72 , 0 ), and 

77eddy = 7? + min(-£i3, -£2^, -£^ 2 ). 

4.2. Numerical results: eddy difFusivity 

Using the same algorithms as employed for R-IV, we 
have computed the eddy diffusivity tensor (see Fig. 3) 
for mTG for a = 6 = 1, as in Lanotte et al. (1999), 
the coefficient c and the molecular diffusivity 77 rang¬ 
ing in the intervals [0.25,0.4] step 0.05 and [0.1,0.16] 
step 0.001, respectively. Advective terms were com¬ 
puted by pseudospectral methods with the resolution 
of 48^ Fourier harmonics. Dealiasing was performed 
by keeping in the solution only harmonics with wave 
numbers not exceeding 21. Energy spectra decaying 
by 7-10 orders of magnitude, this resolution is suffi¬ 
cient. As in the case of R-IV, iterations were termi¬ 
nated, when an estimate of the dominant eigenvalue 
was below 10 “^*^ in absolute value and the norm of 
the discrepancy for the normalised associated eigen¬ 
vector was below 5 • 10“^^. Computation of one eddy 
diffusivity correction tensor £> requires 10-20 minutes 
of a 3.9 MHz Intel Core i7 processor (the code is se¬ 
quential). We have also carried out computations for 
the parameter values 

a = 6=1, c=5/13, d = 2/li (45) 

used by Lanotte et al. (1999). Although our algorithms 
and codes are independent from those applied by Lan¬ 
otte et al. (1999), our values of 77eddy coincide with 
those found ibid, in four significant digits. 

In all runs shown in Fig. 3, we have found that 
£23 > 0 and £§3 < 0, £^3 < 0. Thus, negative eddy 
diffusivity gives rise to growing large-scale magnetic 
modes with horizontal wave vectors of the large-scale 
harmonic modulation. Physically the most interest¬ 
ing case occurs when generation of large-scale fields 
is not obstructed by generation of small-scale fields. 
The segments of the plots of the minimum eddy diffu- 
sivities corresponding to this case are shown by bold 
solid lines. Each segment is bounded on the left by the 
critical point for the onset of generation of small-scale 



Eigure 3. Minimum eddy diffusivity 77eddy = V ~ £23 
(vertical axis) in mTC (39) for a = 6 = 1. Bold 
solid lines: the segments of plots for molecular dif- 
fusivities 77 , for which a large-scale, but no small-scale 
magnetic field is generated. Lower to upper curves: 
c = 0.25, 0.3, 0.35, 5/13, 0.4. 



Eigure 4. Growth rates (vertical axis) of dominant 
small-scale magnetic eigenmodes (solid line: the AAAEO 
subspace; dashed line: the SAADE/ASAOO subspaces) 
generated by mTG (39), (45). 

magnetic field, and on the right by the point where 
eddy diffusivity becomes positive. The critical values 
of molecular diffusivity for the onset of generation of 

Table 1. Critical molecular diffusivities 77 for the onset 
of generation of small-scale magnetic fields in three 
invariant symmetry subspaces. 


c 

AAAEO 

SAAOE/ASADD 

SSAEE 

0.25 

0.1143 

0.1580 

0.1203 

0.3 

0.1091 

0.1491 

0.1118 

0.35 

0.1065 

0.1333 

0.0985 

5/13 

0.1071 

0.1105 

0.0849 

0.4 

0.1080 

0.0851 

0.0770 
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small-scale magnetic fields in some invariant subspaces 
(see the previous section) are shown in Table 1. The 
dominant magnetic eigenmodes have been computed 
applying the algorithms of Zheligovsky (1993) with a 
resolution of 64^ harmonics, the dealiasing was per¬ 
formed by keeping harmonics with wave numbers not 
exceeding 29. Energy spectra of the obtained eigen¬ 
modes decay by at least 11 orders of magnitude. For 
the considered rj, the dominant eigenmodes belong to 
the AAAEO subspace (see Fig. 4; we did not aim at 
computing the dominant magnetic eigenmodes in all 
symmetry subspaces). The dominant eigenmodes in 
the AAAEO and SSAEE subspaces turn out to be y-sym- 
metric. The plots of ? 7 eddy have vertical asymptotes 
located at the critical values for the onset of gener¬ 
ation of the small-scale magnetic field in the SSAEE 
subspace (see Zheligovsky 2011 for explanations). 

4.3. Numerical results: finite scale separation 

We now consider the case of a finite (i.e., non¬ 
infinitesimal) scale separation e. By comparing nu¬ 
merical solutions with the multiscale predictions, we 
can roughly estimate the range of the scale ratios e, 
for which the asymptotic formalism qualitatively cor¬ 
rectly describes the large-scale dynamo driven by an 
array of mTG flow cells. As established in the pre¬ 
vious section, for a high scale separation (i.e., in the 
limit of small e), a large-scale magnetic mode gener¬ 
ated by mTG grows the fastest, when the unit wave 
vector q is horizontal and B = ea in the large-scale 
modulation (18). Such a mode is asymptotically close 
to 

b = Re(e-‘i’^(s 3 (x)+e 3 (46) 

2 

+ is qmGmsi^) + . 

m—1 

To study directly magnetic field generation for an 
arbitrary finite scale separation e, we can employ the 
procedure used by Zheligovsky et al. (2001). Namely, 
we consider the problem ( 8 ) for a field of the form 

b(x) =e‘"‘i'’‘b'(x), (47) 

where q is a constant unit wave vector. A small-scale 
(i.e., having the spatial periodicity cell T^) vector field 
b'(x) satisfies the eigenvalue equation® 

£,qb' = - s^lqpb') + Vx X (v X b') (48) 

-I- ie(q X (v X b') + 2r]{q- Vx)b') = Ab' 

and the corollary of the solenoidality condition 

Vx • b' -I- ieq • b' = 0. (49) 

®We have preserved the factor |qp in (48) for this equation to 
remain valid for any vector q, and not just for a unit one. 


This approach is advantageous in that it does not re¬ 
quire performing the asymptotic analysis of Section 2 
and is applicable for all scale ratios s, and not only very 
small ones. However, it is less general in that, on the 
one hand, a solution to the eigenvalue problem (48)- 
(49) provides information for only one instance of the 
amplitude modulation vector eq. On the other, it is 
only applicable when tackling a linear stability prob¬ 
lem such as the kinematic dynamo problem studied 
here, but does not deliver a simplified statement of a 
weakly nonlinear stability problem. 

For £ > 0, even and odd vector fields (that are lin¬ 
ear combinations of Fourier harmonics such that the 
sum of the wave numbers in the horizontal directions 
is even or odd, respectively) constitute invariant sub¬ 
spaces of iieq (48). If q = for i ^ m, vector fields, 
symmetric or antisymmetric in Xm, also constitute in¬ 
variant subspaces. The case z = to is more subtle: 
vector fields, whose real part is symmetric or antisym¬ 
metric in Xm, and the imaginary part is, respectively, 
antisymmetric or symmetric in Xm, constitute two in¬ 
variant sets. However, these sets are not linear sub¬ 
spaces (over the field of complex numbers); in other 
words, this property can be used in computations, but 
it does not restrict an eigenmode, since multiplying an 
eigenmode by the complex unity i does not give rise 
to a new eigenmode — except for £ = 0 , when only 
the symmetric or antisymmetric part of the eigenmode 
“from which the branch originates” is non-zero. Gon- 
sequently, for q = we can use labels for branches of 
dominant eigenfields of Xleq that have the same mean¬ 
ing as the labels of invariant subspaces of the domain 
of the small-scale magnetic induction operator ex¬ 
cept for the symmetry or antisymmetry in place to 
of the label is determined only for the eigenmode for 
£ = 0. The symmetry 7 , involving swapping of the 
horizontal Gartesian coordinates as well as swapping 
of vector field components, does not distinguish invari¬ 
ant subspaces of iieq for £ > 0. It maps eigenfunctions 
of iieq to eigenfunctions of iigq' for q' = ( 92 , 9i, 93 )- 

We have computed the dominant eigenvalues (i.e., 
the ones having the maximum real part among all 
eigenvalues for the given parameter values) of the mag¬ 
netic induction operator and the associated large-scale 
magnetic modes generated by mTG (39), (45) — the 
flow employed by Lanotte et al. (1999) — for the 
wave vectors of the large-scale amplitude modulation 
q = (1,0,0) (see Figs. 5 and 6) and q = (l,l,0)/-\/2 
(Figs. 7 and 8 ). Since the flow possesses the symme¬ 
tries in xi and X 2 and the 7 -symmetry, actually the 
computations cover all possible choices of q from the 
following list: ±ei, ± 02 , (± 1 , ± 1 , 0 )/-\/ 2 . 

Plots of growth rates of large-scale magnetic modes 
for q = ei and a varying scale ratio £ are shown in 
Fig. 5 for ?7 = 0.1 used by Lanotte et al. (1999), as 
well as for rj = 0.11 and 0.12. For these molecular 
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Figure 5. Growth rates (vertical axis) of dominant 
large-scale magnetic modes (SSAEE subspace) gener¬ 
ated by mTG (39), (45) for q = ei. 



Figure 6. Growth rates (vertical axis) of domi¬ 
nant large-scale magnetic eigenmodes generated by 
mTG (39), (45) for rj = 0.02, q = ei: four branches 
in symmetry subspaces AAAED (bold line), AAAEE (thin 
dotted line), SSAED (bold dotted line and thin solid 
line for 0.73 < e < 0.78: the real and imaginary parts 
of the associated eigenvalues; outside this interval, the 
eigenvalue is real), and SSAEE (thin solid and dashed 
lines: the real and imaginary parts of the associated 
eigenvalues). 

diffusivities the dominant eigenvalues of the operator 
£eq are real. Zheligovsky et al. (2001) noticed that 
a graph of the dominant growth rates is periodic in e 
with period 1 (because any large-scale field b'(x), 
where b'(x) is a small-scale field, can be also expressed 
as (e*^'i ’^b'(x)), and for an arbitrary inte¬ 

ger p the field b'(x) is also small-scale). Also, a 
graph of the dominant magnetic mode growth rate as 
a function of the scale ratio e is symmetric about the 
vertical axis: applying complex conjugation to equa¬ 
tions (48) and (49) shows that if, for a given scale ratio 
e, b'(x) is a small-scale eigenfunction associated with 



Figure 7. Growth rates (vertical axis) of dominant 
large-scale magnetic modes generated by mTG (39), 
(45) (solid lines) for q = (1,1,0)/-\/2. For comparison, 
growth rates of dominant large-scale magnetic modes 
for q = ei and the same molecular diffusivities rj are 
shown (dotted lines; cf. Fig. 5). 



Figure 8. Growth rates (vertical axis; solid line) of 
dominant large-scale magnetic eigenmodes generated 
by mTG for rj = 0.02, q = (1,1,0)/-\/2. For compar¬ 
ison, growth rates in the branches of dominant large- 
scale magnetic modes for q = ei in the symmetry 
subspaces AAAED (right) and SSAEE (left) for the same 
rj are shown (dotted lines; cf. Fig. 6). 


an eigenvalue A, then b'(x) and A are, respectively, 
a small-scale eigenfunction and the associated eigen¬ 
value for the opposite ratio —e. Consequently, graphs 
of the dominant growth rate are symmetric about each 
vertical line s = qj2 for integer q. By contrast, the 
plots in Figs. 5 and 6 show eigenvalues associated with 
branches of eigenfunctions of the problem (48)-(49), 
smoothly parameterised by e. They have a period 2 
in e and are symmetric about vertical lines s = q for 
all integer q. The parabolic shape of the plots near 
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£ = 0 agrees with expansion (9) for Aq = Ai = 0. That 
£ = 0 is a local minimum of the plots in Fig. 5 cor¬ 
roborates that magnetic eddy diffusivity is negative 
for the molecular diffusivities rj = 0.1, 0.11.0.12, for 
which plots are presented in this figure; the respective 
eigenmodes b'(x) constitute SSAEE branches. 

Near the origin, the plots of growth rates in Figs. 5 
and 7 have a parabolic shape (which is a signature of 
magnetic eddy diffusivity) for e below 0.1; this roughly 
estimates the range of scale ratios for which the asymp¬ 
totic formalism describes qualitatively correctly the 
large-scale dynamo driven by an array of mTG flow 
cells. A similar parabolic-shape correction of growth 
rates due to the action of eddy diffusivity is observed 
for non-neutral magnetic modes (Figs. 6 and 8) in all 
the symmetry subspaces considered. 

Our computations demonstrate that for rj = 0.1, 
mTG can generate large-scale magnetic field by the 
mechanism of negative eddy diffusivity in a range of 
parameter values. By contrast, for 77 = 0.02 no large- 
scale magnetic field generation was found by Devlen et 
al. (2013) in DNS. We have computed four branches 
of dominant eigenmodes for 77 = 0.02 and q = ei 
(see Fig. 6), that belong to invariant subspaces AAAED, 
AAAEE, SSAEO and SSAEE with the resolution of 96^ har¬ 
monics (upon dealiasing, harmonics with wave num¬ 
bers up to 45 are kept); energy spectra of the eigen¬ 
modes decay by at least 9 orders of magnitude. 

We observe two major differences with the case 
77 = 0.1. First, a small-scale dynamo persists for 
77 = 0.02. Implementation of the TFM procedure 
requires integrating equation (2); the solution con¬ 
verges to the dominant small-scale mode, amplitude- 
modulated by the large-scale harmonic e"^^™. Glearly, 
in the presence of a small-scale dynamo, the solution is 
dominated by the growing small-scale mode, and not 
by the neutral mode (46). Solutions can be expanded 
in the series (9) in the scale ratio e, the series for the 
eigenvalue A now beginning with the respective small- 
scale dynamo eigenvalue. For a parity-invariant flow 
this modihes the molecular diffusivity operator, acting 
on the amplitude-modulating factor (called amplitude) 
in the respective large-scale mode; like in the absence 
of a small-scale dynamo, the correction is due to in¬ 
teraction of the fluctuating part of the magnetic field 
and the small-scale flow, and thus again eddy diffusiv¬ 
ity is the leading-order eddy effect. Second, the point 
£ = 0 is now a local maximum, implying that eddy dif¬ 
fusivity is now positive. However, the growth rates of 
large-scale magnetic modes are still positive when |e| 
is small, i.e., these modes do grow, albeit slower than 
the small-scale modes for e: = 0. In other words, the 
growing large-scale modes decay relative to the faster 
growing small-scale mode, which explains the state¬ 
ment “a dynamo is observed but it is not a large-scale 
dynamo” (Devlen et al. 2013). 


Yet another difference with the case 77 = 0.1 is vis¬ 
ible in the behaviour of dominant eigenmodes consti¬ 
tuting the SSAEE branch. For 77 = 0.02, they expe¬ 
rience two bifurcations: on increasing the scale ratio 
s, a pair of real eigenvalues (including the dominant 
one) turns into a pair of complex-conjugate ones at 
e « 0.34, that are superseded again by two real eigen¬ 
values at e « 0.98 (only the largest of which is shown in 
Fig. 6). We observe a characteristic feature of depen¬ 
dence on the parameter near a point of such a bifurca¬ 
tion: the plots of real eigenvalues and of the imaginary 
part of complex eigenvalues (but not of the real part 
of the complex eigenvalues) have singularities of the 
kind of ^/x for x > 0 near zero — the growth rate de¬ 
pends on £ continuously, but its derivative is infinite. 
This stems from the fact that the quadratic charac¬ 
teristic polynomial of Xleq, reduced onto the invariant 
plane of the associated eigenfunctions, has coefficients 
that are differentiable in £, and hence the discriminant 
is approximately a linear function of £ near the point 
of bifurcation. Vanishing of the discriminant at such 
points gives rise to the singularities mentioned above 
(the real parts of complex eigenvalues are not affected, 
since they are just proportional to the coefficient of the 
linear term of the characteristic polynomial). Gher- 
tovskih et al. (2010) observed a similar behaviour in 
the dependence of magnetic field generation by ther¬ 
mal convection on the rotation rate (see Fig. 18 ibid.). 

We have also computed the short-scale parts b'(x) 
of the dominant large-scale magnetic modes (47), gen¬ 
erated by the same instance of mTG (39), (45) for the 
wave vector q = (l,l,0)/-\/2 and molecular diffusiv¬ 
ities 77 = 0.1, 0.11, 0.12 (using the resolution of 64^ 
Fourier harmonics) and 0.02 (96^ harmonics). As for 
q = ei, this has been done by solving the eigenvalue 
problem (48)-(49) for the modified operator of mag¬ 
netic induction Xleq. For all considered 77 and £, the 
computed dominant short-scale modes of iieq possess 
now the 7-symmetry, the antisymmetry in X3 and the 
symmetry about the xa-axis, which is the composition 
of the symmetries in xi and X2: 

f{-xi,-X 2 , X 3 ) = -/Xx), 1 = 1,2, 

fi-Xi,-X2,X3) = /^(x). 

These short-scale modes are comprised of the Fourier 
harmonics, for which all the three wave numbers ki in 
the directions Xi have the same parity. The associated 
eigenvalues of the operator Xf^q are real. 

For mTG, eddy diffusivity is the same for all hor¬ 
izontal wave vectors (see (43)). Comparison of the 
eigenvalues computed for q = (1,1,0)/v^ and q = ei 
in Figs. 7 and 8 illustrates how this axisymmetry is 
reflected in the eigenvalues for £ > 0. We observe that 
the dependence of the dominant eigenvalues on the di¬ 
rection of a horizontal wave vector is very weak when 
£ is as large as roughly 0.8 for 77 = 0.1, when £ < 0.7 
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for rj = 0.11 and 0.12, and only when e < 0.22 for 
r] = 0.02. 

5. TFM vs MST: analytic and numerical 
comparison 

We have seen in Section 3 that the TFM used by 
Devlen et al. (2013) for evaluating magnetic eddy dif- 
fusivity for R-IV yielded the results compatible with 
those obtained by employing the homogenisation tech¬ 
niques within the MST approach. Given that distinct 
types of averaging are employed in MST and TFM, 
this conformity of results may seem unexpected. In 
the present section we compare the two approaches. 

TFM starts by computing a zero-mean solution b' 
to equation (2) for the test field 

btest = e‘'^®’"e„ (50) 

(this is equivalent to employing the two real fields 
(6) for k = EBjn, but simplifies the algebra). Any 
solenoidal small-scale zero-mean field (for instance, 0) 
can serve as an initial condition for b'. The solution 
will then automatically be solenoidal at any time t > 0. 
TFM assumes that (2) does not have growing solutions 
for the test fields and averaging applied. Numerical in¬ 
tegration of (2) proceeds till transients decay and the 
solution b' saturates. The eddy diffusivity correction 
tensor is then deduced as the matrix that relates the 
obtained mean e.f.m. v' x b' with the test fields (50). 

5.1. TFM with volume averaging 

We now consider a variant of TFM, in which volume 
averaging is involved in extracting the fluctuating part 
of the auxiliary fields, b', and show that then the TFM 
values of eddy quantities converge in the limit of large 
scale separation to the values yielded by MST. The 
demonstration, given here for steady flows v', can be 
readily extended to encompass time-periodic flows. 

Our solutions can be obtained as the real and imag¬ 
inary parts of the fluctuating part of the field 


£ > 0, the elliptic operator Xige™ the r.h.s. of (52) 
is an 0(e) perturbation of the operator of magnetic 
induction, XI. Consequently, this stage of TFM can 
be readily understood in the framework of MST. By 
the general theory of perturbation of linear operators 
(Kato 1966; see also Vishik 1987), an eigenfunction of 
XI and the associated eigenvalue involved in a Jordan 
cell of size M are altered by 0(e^/'^). 

TFM is applicable when no small-scale dynamo op¬ 
erates. In this section we assume that the kernel of 
the operator of magnetic induction, defined in the box 
of periodicity of the flow, is three-dimensional (for a 
given V, this holds for all r/ > 0 except only for a count¬ 
able number of rj values). A solution to (52) is a sum 
of a transient b*'', whose rate of exponential decay is 
0(1), and the neutral mode of the perturbed operator, 
that branches from the respective neutral mode of XI 
(for which M = 1): 

b = S„(x) -I- e„ -I- 0(e) -I- b**' (53) 

for any permissible initial conditions for b. 

1. Magnetic a-effect. For a generic steady flow v, 
we can now calculate the TFM estimate of the a-tensor 
using the ansatz (4). By (51) and (53), after the tran¬ 
sient decays below 0(e) at times O(lne), 

b = e*^“'™(S„(x)-I-e„-I-0(e)). (54) 

Large-scale computations of b' are usually done for 
a rational e = ii/i 2 (with common factors cancelled 
out in integers zi and 12 ) such that the periodicity of 
gieajm jg compatible with that of the small-scale flow v. 
Thus, we can assume that the computational domain 
has the size 27rz2 in Xm- When applied to a steady field, 
the Fourier transform (5) involved in (4) differs only 
by a constant factor from the inverse Fourier transform 
= (l/f^)-^k,o that recovers coefficients in expansion 
of a function in the spatial variables: 

=/k- 


b = btest + b' = e'^'^^b, (b) = e„. (51) 

Note that when extracting the fluctuating part b', we 
average b after pulling out the factor e*^^™, since aver¬ 
aging over Xm any field of the form f (x) , where f 
is independent of £, yields just 0. The evolution equa¬ 
tion (2) for the auxiliary field b' is then equivalent to 
the equation obtained by substituting (51) into (2) and 
cancelling out the exponential: 

^ = (-Cee„b)', (52) 

(the operator Xiee„ is defined by (48)); b also satis¬ 
fies the condition (49), stemming from solenoidality of 
b', and has a constant average (b) = e„. For small 


Here V denotes the volume of the spatial periodicity 
domain. Using (54), we find 


1 * 277^2 


/*27r /*27r 


(v' X b') = 


X b' 


=(v' X S„) + 0(£) 


dx 

(27r)3z2 

(55) 


(here any spatial averaging is acceptable, provided 
it does not involve averaging in Xm, or otherwise spe¬ 
cial precautions are taken as discussed above). Since 
(btest) = e„, by (4) the n-th column of the 3x3 
matrix a coincides in the limit e —> 0 with the n-th 
column of 21 . 

Remark 1. TFM for evaluation of the a-effect ten¬ 
sor in non-parity-invariant flow in the original formu¬ 
lation (Schrinner et al. 2005, 2007) prescribed the use 
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of constant test fields btest = which coincides with 
(50) for £ = 0. Consequently, iige™ = -C in (52), and 
thus this version of TFM with the spatial averaging 
reproduces the MST a-effect tensor precisely. 

2. Magnetic eddy diffusivity. If the flow is parity- 
invariant, i.e., v(—x) = —v(x), the three small- 
scale eigenfunctions from the kernel of £ are parity- 
antiinvariant: S„(—a;) = S„(a:). The parity-invariant 
part of b, even if zero initially, is subsequently pro¬ 
duced from the predominantly parity-antiinvariant 
held (53) by the term iee™ x (v x b) in (52). Since all 
parity-invariant eigenmodes of S, decay (by the orig¬ 
inal assumption on the spectrum of £,), the parity- 
invariant part of b remains 0(e) at all large enough 
times. We can seek b as a perturbed truncated series 
for the neutral mode of £, known from MST: 

b = e„ + S„ + ieG™„ + b^"'. 

Substituting this ansatz into (52), we obtain an equa¬ 
tion of the form 


small-scale held S„(x) and at least one period of the 
modulating harmonic 

Above, we have investigated the algorithms for eval¬ 
uating for n ^ m. Test helds (50) for m = n are 
gradients and hence incompatible with our analysis. 
To evaluate 2D„„, we can use solutions to (2) for test 
helds, that are real and imaginary parts of 

btest = (ie)-'V X 

where ni ^ n, n 2 ^ n and j is an integer. 

5.2. TFM with other spatial averagings 

We now consider briehy the canonical variants of 
TFM, in which the averaging denoted by a bar is per¬ 
formed over one or two Cartesian variables under the 
same assumptions as in Section 5.1. For simplicity, 
we again ignore the memory effect by assuming that 
the test held does not depend on time. As before, we 
cancel out in (2) the exponent involved in the 

unknown held (51), and hnd 


^bnew 

dt 




0(£^). 


We therefore hnd 


b — e‘®“’"(e„ -|- S„ -|- ieGmn + 0(£^) + b**^), (56) 

where b*’’ is a transient, whose rate of exponential de¬ 
cay is 0(1). 

We now calculate the entries of the magnetic eddy 
diffusivity correction tensor from the equation 

S'£e„(v' X b') = (57) 

p,q ^ P ''q 


This is ansatz (4) for steady how and zero a-effect. 
As in item 1 , we assume e = ii/z2 is rational so that 
the periodicities of b (56) and the small-scale how v 
are compatible, and the computational domain has the 
size 27ri2 in Xm- On the one hand, we then hnd 


(v' X b') 



(v' X (S„ -k ieG mn )) + 0{e 


dx 

(27r)3i2 


k2)„„ + 0(e2). (58) 


On the other. 


By (57), iDmn = -lime^oi7mn for all m ^ n, i.e., 
TFM does produce in the limit £ —> 0 the respective 
entry of the tensor of magnetic eddy correction. Note, 
however, that a sufficiently high spatial resolution is 
necessary for a satisfactory discretisation of both the 


— =ipi:ee„b. (59) 

Here denotes a projection that deletes the mean 
held, but preserves the volume average: 

qjf = f-f-k (f). 

While (2) is equivalent to (59), the latter equation 
has advantages: (z) It can be numerically integrated 
in the how periodicity cell without encountering 
the instabilities of problem (59) which may exist at 
larger spatial scales (note that computations must be 
done in a box of size 27r/£ in Xm when the exponential 
or sinusoidal dependence on exm is preserved in b'). 
Such instability will then manifest itself by unbounded 
amplihcation of the growing eigenfunctions of the op¬ 
erator that emerge from round-off errors, and 

this will progressively wipe out the contribution from 
the inhomogeneity in (2) which we are looking for. 
(zz) It enables us to compute auxiliary helds b' for 
irrational £ without suhering from problems due to 
the presence of two incommensurate spatial frequen¬ 
cies in the solution, (zzz) One can apply to solutions of 
(59) spatial averaging over any variable, including Xm- 
In turbulence computations, which are made in the 
large (from the prospective of the present discussion) 
computational box, averaging a held after cancelling 
out the exponential is also a feasible operation that is 
just equivalent to computing the appropriate Fourier 
transform. 

For any permissible initial conditions, (59) admits 
solutions similar to (53): 

b = e„ -|- s„ -|- 0(£) + b*"*^ 
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and, for parity-invariant flows, similar to (56): 

b = e„ -h s„ -h iegmn + 0 (£^) -K b*"', 

where b**^ are transients, whose rate of exponential 
decay is 0(1). The fields s„ and gmn have zero av¬ 
erages: s„ = = 0 , the fields e„ -|- s„ belong to 

the kernel of the operator For parity-invariant 

flows, s„ are parity-antiinvariant and g^n are parity- 
invariant. But here the similarity ends, e.g., s„ ^ S(j 
and gmn 7 ^ Consequently, in the limit e —>■ 0 we 

can expect a qualitative but not quantitative agree¬ 
ment of MST results with those of TFM with a non¬ 
volume averaging. 

Remark 2. Plane-parallel flows independent of a 
Cartesian coordinate Xm are a special case, for which 
a solution (51), harmonically modulated by the factor 
involves the small-scale part b that is indepen¬ 
dent of Xm- Consequently, for such flows, s„ = S„ and 
gjn = Gjn for j ^ m ^ n, and hence TFM recovers 
precisely the components 'Z>jn of the eddy correction 
tensor. This is the case of R-IV. 

5.3. Kinematic generation by mTG: magnetic 
structures 

To understand the absence of negative magnetic 
eddy diffusivity in the TFM results of Devlen et al. 
(2013), we first inspect magnetic modes obtained from 
numerical solutions of the underlying eigenvalue prob¬ 
lem for mTG (39), (45) and 77 = 0.1. The modes are 
eigenfunctions of the magnetic induction operator and 
give rise to exponential in time solutions of the mag¬ 
netic induction equation 

9b 

- = i:b. (60) 

We consider first magnetic eigenmodes with the peri¬ 
odicity box of size (27r)^. As discussed in Section 4.1, 
due to the symmetries of the flow, magnetic modes 
have symmetries or antisymmetries in Cartesian coor¬ 
dinates Xi and in each mode the sums of wave numbers 
ki -f k 2 in all harmonics have the same parity, as well 
as all sums fci -I- fca. These 5 symmetries are indepen¬ 
dent and split the domain of the magnetic induction 
operator into 32 invariant subspaces. On top of this, 
magnetic modes can be symmetric or antisymmetric 
with respect to swapping of the horizontal coordinates 
xi -fT- X 2 (the symmetry 7 ), but this symmetry is not 
independent of the 5 former ones: it splits into invari¬ 
ant subspaces only 8 of the 32 aforementioned invari¬ 
ant subspaces — namely those, in which the sums of 
wave vectors fci -I- k 2 are even, and vector fields are ei¬ 
ther symmetric in both xi and X 2 ^ or antisymmetric in 
both of these Cartesian variables. Thus, the symme¬ 
tries of mTG split the domain of the magnetic induc¬ 
tion operator into 40 invariant subspaces. We have 


computed dominant (i.e., having the largest growth 
rates) magnetic modes in each of them. 

Only in 3 subspaces out of 40, growing 27r-periodic 
magnetic modes have been found (see Table 2). We 
first inspect suitably averaged fields; as discussed in 
the next section a particularly revealing average is that 
over the X 3 coordinate, the mean field being a func¬ 
tion of Xi and X 2 - In Fig. 9 we show such mean fields 
63 ( 0 : 1 , X 2 ); clearly, they do not survive horizontal aver¬ 
aging over the {xi,X 2 ) plane, because the positive and 
negative contributions in Fig. 9 cancel. This figure also 
illustrates some of the symmetries of the invariant sub¬ 
spaces, to which the 3 modes belong. Figures 10 and 
11 show isosurfaces of the energy at the level |bp = 2 , 
and of the vertical magnetic component at the level 
63 = 2/3, for two dominant modes, that are not mutu¬ 
ally related by any of the symmetries. (The dominant 
modes in subspaces SAAOE and ASADO are mapped into 
each other by the symmetry 7 , and the mode in AAAEO 
is 7 -symmetric.) 

In slow dynamos, magnetic structures can be re¬ 
lated to stagnation points of the flow. Eight families of 
stagnation points of mTG are listed in Table 3; we have 
checked numerically that no other stagnation points 
exist in mTG (39), (45). Each of the first four families 
is a 7 -symmetric set; families V and VI are mapped 
by 7 into each other, as well as families VII and VIII. 
Lines joining stagnation points of family I and paral¬ 
lel to Cartesian axes constitute a heteroclinic network: 
any such vertical line consists of heteroclinic trajecto¬ 
ries connecting adjacent stagnation points of families 
I and II, and a horizontal line consists of heteroclinic 
trajectories connecting a pair of adjacent stagnation 
points of family I. Each plane, parallel to a Cartesian 
coordinate plane and containing stagnation points of 
family I, is cut by the aforementioned heteroclinic tra¬ 
jectories into squares of size tt, which are invariant sets 
for mTG (this stems from the proportionality of Vi to 
sinxi for each i). Vertical and horizontal lines join¬ 
ing stagnation points of family IV constitute another 


Table 2. Maximum growth rates. A, of 27r- and 47 r- 
periodic magnetic modes with different symmetries 
generated by mTG (39), (45) for ry = 0.1. 


Period 

Symmetry subspace 

A 

27T 

SAAOE,ASADO 

0.01602 

27T 

AAAEO 

0.01383 

4:TT 

AAAOOE,ASAOOE, SAADOE, SSADDE 

0.01763 

4:TT 

ASAOEE,SAAEOE,SSAEOE, SSADEE 

0.01734 

4:TT 

ASAEEE,SAAEEE 

0.01602 

4:TT 

AAAOEE,AAAEOE,ASAEOE, SAADEE 

0.01404 

4:TT 

AAAEEE 

0.01383 

4:TT 

AAADDO,SAAOOD,ASAODD,SSAOOO, 
AASOOD,SASOOO,ASSDOO, SSSDDO 

0.00226 
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Figure 9. The xa-averaged normalised mean fields b^{x\,X 2 )/for 77 = 0.1 in the box periodicity of size (27r)^. 
The same colour-coding scheme is used in all panels; the data outside the interval [—5, 5] is clipped. 



Figure 10. Isosurfaces of the energy for two rms-normalised dominant 27r-periodic modes: SAAOE, |bp = 2 (left), 
AAAEO, |bp = 2 (centre) and AAAEO, |bp = 2/3 (right). 



Figure 11. Isosurfaces of the vertical magnetic field component, 63 = 2/3, for two rms-normalised dominant 27r- 
periodic modes, ASAOO (left) and AAAEO (right). 


18 





















Table 3. Stagnation points of modified Taylor-Green flow (39) for a = 5 = 1, and the spectral structure (eigenvalues, 
(7, and the proper subspaces) of the Jacobian matrix, [dv^/dxn], at these points; ji and j are arbitrary integers. 


Family 

Stagnation point 

Proper subspace 

Eigenvalue, a 



ei 

2 -k(-l)^'i+^'"+^ 3 ( 2 -k 3 c) 

I 

7 !'(jl,j 2 , js) 

©2 

2 - (- 1 ) 31 + 32 + 33(2 -p 3 c) 



63 

-4 



ei 

-2 

II 

+ 5) 

e2 

-2 



63 

4 

III 

’’■(jl + 5 ,J 2 + 5,73) 

{61,62} 

—2 ± i \/c -k 2 



63 

4 



61 

2 

IV 

7r(ji + hh + 5V3 + 5) 

62 

2 



63 

-4 

V 

’’■(jl + 5,J2 ,43 + 5) 

crei - 4 (-l)Ji+J 2 +*( 3 c- l)e3 

l±v/l + 4 ( 3 c-l)( 2 -c) 



62 

-2 

VI 

+ 5,43 + 

crea -k 4 (- 1 )Ji+J 2 +J 3 ( 3 c _ 1)03 

l±A/l-k 4 ( 3 c-l)( 2 -c) 



61 

-2 

VII 

(a:i,j7r,a:3), where 

61(0- - C 5 (a:i)) -k 63(0(311,4:3) 

(3(311,313) ± VC 4 ( 31 l, 3 l 3 ) 


sin^ a;i = C 2 "^ (“Ci ± VCi + 4C2) 

62 

-2(3(311, X3) 


cosxs = (—l)^( 3 c— 1 ) sinxi tanxi 



VIII 

(j 7 r,X 2 ,X 3 ), where 

61 

-2(3(312, X3) 


sin^ X 2 = (2 (-Cl ± \/Ci + 4C2 j 

62(0- - C 5 (a: 2 )) + 63(0(312,3:3) 

(3(312,33) ± 1/(4(32,33) 


cos 0:3 = (— 1 )^ (1 — 3 c) sina;2 tan 0:2 




Note. Here G = 3c(3c + 1), <^2 = 8(3c — l)(2c — 1), ^ 3 ( 0 ;, z) = (—l)t((3c — l)((3c/2 + l)sin^x — 2sin^a;) — cos 22 ), 
Ciix) = ((—l)'^(6c — 2) — 4c + ((—l)t(3c — 27c^) + 4 — 14c + 54c^) sin^ a:)/(2c — 1), (q{x,z) = 4(3c — 1)(2 — sin^ a;) sin a; sin 2 , 
C,i{x, 2 ) = (2 sin^ 2 — (8c + (5 — 19c + 3c^ + 90c^) sin^ x)/{l — 2c))^ + 4(1 — 3c)(2 — sin^ a;) sin^ a;sin^ 2(2 + 3c — (8 — 20c) sin^ x) for 
odd j and Q 4 ,{x, 2 ) = (2cos^ 2 + (—3 + 3c + 9c^) sin^ a:)^+ 4(1 — 3c)(2 — sin^ x) sin^ a;sin^ 2(2 + 3c — (8 — 20c) sin^ a:) for even j. 


heteroclinic network: they consist of heteroclinic tra¬ 
jectories connecting points of family IV with adjacent 
stagnation points of families III, V and VI. 

The Jacobian matrix of a solenoidal flow generically 
has either one positive eigenvalue and two eigenval¬ 
ues with negative real parts, or one negative eigen¬ 
value and two eigenvalues with positive real parts. In 
the vicinity of a stagnation point of the former kind 
(having a one-dimensional unstable manifold), mag¬ 
netic flux ropes usually emerge (Moffatt 1978; Gal¬ 
loway & Zheligovsky 1994) that are aligned with the 
unstable direction. Near a stagnation point of the lat¬ 
ter kind (possessing a two-dimensional unstable man¬ 
ifold), magnetic sheets typically emerge (Ghildress & 
Soward 1985) spreading along the unstable manifold. 
(Formation of these magnetic structures may be pro¬ 
hibited by symmetries.) 


We observe such patterns of asymptotic nature, 
foremost, vertically oriented flux ropes, that are cen¬ 
tred at stagnation points of family III (whose one¬ 
dimensional unstable manifolds are segments of ver¬ 
tical lines), in the plots of isosurfaces of the magnetic 
energy at the level |bp = 2 (Fig. 10, left and cen¬ 
tral panels) and of the vertical component of magnetic 
field (Fig. 11) for both modes, shown in the figures, 
from the symmetry subspaces SAAOE and AAAEO. These 
“principal” ropes terminate near stagnation points of 
family IV, whose two-dimensional unstable manifolds 
are horizontal planes, and which give rise to magnetic 
field sheets revealed by energy isosurfaces at the low 
level |bp = 2/3 (Fig. 10, the right panel). The sheets 
intermix into vertical flux ropes centred at stagnation 
points of family II. In the AAAEO mode, adjacent prin¬ 
cipal flux ropes are oppositely directed (see Fig. 9); 
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consequently, the flux ropes associated with stagna¬ 
tion points of family II are comprised of two pairs 
of oppositely oriented “flux fibres” (such compound 
flux ropes were considered by Galloway & Zheligovsky 
1994). Since fine structures are accompanied by en¬ 
hanced dissipation, the compound ropes are weak and 
not seen in the right panel of Fig. 10 — these rela¬ 
tively high-level isosurfaces only determine the region 
in space, where the four-fibre flux ropes are located. 
Compound flux ropes consisting of two oppositely di¬ 
rected fibres centred at stagnation points of families V 
and VI are present in the AAAEO mode (these individ¬ 
ual fibres actually look more like beans in the central 
panel of Fig. 10: the width of flux ropes is of the order 
oiRm' , where Rm is the magnetic Reynolds number 
which is clearly not high for 77 = 0.1 considered here, 
and hence the magnetic flux ropes and sheets that we 
observe are rather “fat”). In the SAADE mode, flux 
ropes centred at stagnation points of family VI (but 
not V) are allowed by the symmetries defining the sub¬ 
space; these flux ropes do not have a fibre structure (in 
the left panel of Fig. 10 they are cut into halves by the 
faces of the shown cube of periodicity) and their en¬ 
ergy content is even higher than that of the principal 
ropes. 

All 47 r-periodic magnetic modes, growing for 77 = 0.1, 
are also listed in Table 2. Such modes can be symmet¬ 
ric or antisymmetric in each Cartesian variable xf, this 
is coded by the first 3 characters in the labels (letters S 
and A, respectively) of invariant subspaces, like in the 
case of 27r-periodic modes. The trailing 3 characters 
of the 6 -character labels have now a new meaning: for 
any fixed i, the wave numbers ki in all Fourier har¬ 
monics comprising a 47 r-periodic mode have the 

same parity, which is indicated by letters E or 0 (even 
and odd values, respectively) in position i-|-3. We have 
considered neither the more subtle parity symmetries, 
nor the 7 -symmetry. Since the 6 aforementioned sym¬ 
metries are independent, they split the domain of the 
magnetic induction operator into 64 invariant sub¬ 
spaces. We have computed dominant magnetic modes 
in each of them using 128^ Fourier harmonics (before 
dealiasing), which effectively provide the same spatial 
resolution as 64^ harmonics in computations of the 
27r-periodic modes. 

For the dominant growing 47 r-periodic magnetic 
modes, we show the same plots as for the 27r-periodic 
ones: the mean fields b 3 {xi,X 2 ) averaged over X 3 for 
15 dominant 47 r-periodic modes (Fig. 12), and isosur¬ 
faces of the energy |bp = 2 and of the component 
63 = 2/3 for six of them, that are not mutually re¬ 
lated by any symmetry (Figs. 13 and 14). Clearly, the 
averages over the (xi,X 2 ) plane of the vertical compo¬ 
nent 63 for all dominant modes shown in Fig. 12 are 
zero, as this was the case for the 27r-periodic modes. 
(For the 8 growing modes comprising the last group in 


Table 2, 63 = 0, because they involve only odd wave 
numbers k^.) The most prominent features in Fig. 12 
are the averages of the vertically oriented flux ropes 
centred at stagnation points of family III; all other flux 
ropes cancel out upon averaging over X 3 either mostly 
or completely. It is natural that these mean flux ropes 
of a similar genesis have a similar shape in all panels in 
Fig. 12 and, for instance, have close extremum values, 
the maxima ranging from 5.72 for dominant modes 
from the second group of 47 r-periodic modes in Table 2 
(including subspace ASAOEE) to 6.62 in the fifth group 
(subspace AAAEEE). (The maxima are computed for 
the normalised averages b 3 {xi, X 2 )/■) It turns 
out that the dominant 47 r-periodic modes in subspaces 
ASAEEE, SAAEEE and AAAEEE are just the tiling of the 
cube of periodicity of size Att by 8 cubes of periodicity 
of size 27r with 27r-periodic modes in subspaces ASAOO, 
SAAOE and AAAEO, respectively (note that the growth 
rates of the respective Att- and 27r-periodic modes co¬ 
incide). In fact, each group of 47 r-periodic modes that 
have the same growth rate (see Table 2) are related by 
symmetries. (For instance, the eight slowest-growing 
modes constituting the last group in Table 2 are mu¬ 
tually related by combinations of shifts by 27r along 
the Cartesian axes.) 

5.4. DNS and TFM results for eddy diffusivity 
in mTG 

As noted above, horizontal averaging over the 
(xi,X 2 ) plane cannot be applied to describe a growing 
mean field generated by mTG. Indeed, averaging the 
solenoidality condition for b we find that 63 is spa¬ 
tially uniform at all times; then the spatial average of 
the third component of ( 1 ) shows that it is also time- 
independent. Since 63 cannot grow or decay, such an 
average is unsuitable for studying the negative eddy 
diffusivity dynamo for mTG (for which we are advised 
by MST that 63 7 ^ 0, see (46)). By contrast, planar 
averages can describe growing solutions in the super¬ 
critical case, if one averages along x^ and a diagonal 
direction, or uses any of the two other planar averages, 
over (xijXs) or {x 2 ,X 3 ). Note that, for a flow with a 
large group of symmetries, any planar averaging may 
yield, due to cancellation, identically zero averages for 
modes in certain symmetry subspaces. For instance, 
for mTG, no cancellation occurs for the (a:i,a; 3 ) or 
ix 2 ,X 3 ) averagings for the second and third groups of 
47 r-periodic modes in Table 2, and for diagonal ones 
for the first and fifth groups. Thus, the average over 
(a:i,a; 3 ) or ( 0 : 2 , 3 : 3 ) is adequate in 6 subspaces; the di¬ 
agonal average in 5 subspaces, including the dominant 
one; and none in the remaining 12 subspaces contain¬ 
ing growing modes. None of the easily implementable 
planar averagings is universally applicable. 

We now consider averaging over the ( 372 , 3 : 3 ) plane. 
The evolution of the auxiliary fluctuating field is now 
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Figure 12. The a; 3 -averaged normalised mean field b^{xi, X 2 )/for 77 = 0.1 in a domain of size (47r)^. The white 
dashed lines mark subdomains of size (27r)^. Same colour-coding scheme as in Fig. 9. 
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Figure 13. Isosurfaces of the energy for six rms-normalised isotypic dominant 47 r-periodic modes at the level |bp = 2. 


controlled by the operator ipXlip, where ^ is the pro¬ 
jection that deletes the mean field, but preserves the 
volume average, and £1 is the operator of magnetic in¬ 
duction. The dominant modes of the new operator 
belong to a symmetry subspace, different from those, 
where the dominant modes of XI acting alone reside 
(see in the left panel of Fig. 15 the mean saturated 
magnetic field produced by DNS with the use of the 
Pencil Code^). Despite the additional projections, 
the main visible magnetic structures are still the ver¬ 
tical flux ropes centred at the family Ill stagnation 
points of mTG. The right panel of Fig. 15 shows the 
mean field (which is now a function of xi). It has 
positive and negative extrema at xi = ^71/2. Our 
computations also reveal that the two possible mean 
fields, & 3 (xi) and bz{x 2 ), have the same shape. Again, 
the mean field is anharmonic and therefore the eddy 
diffusivity cannot be spatially constant. 

Owing to the anharmonic nature of the resulting 
mean fields, we must consider test fields involving 


^http://github.com/pencil-code 


many Fourier harmonics. Let us begin with the most 
important contribution from ki = 1. We use again the 
Pencil Code, where TFM is readily implemented. In 
all the cases presented below we have used 72^ mesh 
points. In Fig. 16 we show the results for r] 22 {xi) and 
1733 ( 2 ^ 1 ) for ?7 = 0.1 and fci = 1. Note that both 7711 
and 7733 show strong spatial variations. However, while 
7733 is always positive, 7711 has extended regions where 
it is negative, giving rise to growth of 63 (a;i). 

In principle, negative diffusivities can be used in a 
numerical mean-field simulation. However, one would 
then need to include contributions from larger wave 
numbers ki (or e), where 7722 eventually becomes pos¬ 
itive for large wave numbers. This was demonstrated 
in Devlen et al. (2013), where the turbulent diffusiv¬ 
ity kernel was spatially constant, and so the relevant 
eigenvalue problem became 

AA2 = -(77-I-7722(^:1, iA))A:iA2 

(cf. (16)). Here, A 2 is the Fourier amplitude and, 
for consistency (cf. (4)-(5)), the eddy correction 7722 
should be calculated for w = iA, which is in gen- 
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Figure 14. Isosurfaces of the vertical magnetic component for six rms-normalised isotypic dominant 47 r-periodic modes 
at the level 63 = 2/3. 


eral complex. In the present case, we only find non- 
oscillatory growth, so A is real and therefore the fre¬ 
quency u), for which 7722 is needed, is purely imagi¬ 
nary. Since the dependence of 7722 on ui is in general 
nonlinear, one has a nonlinear eigenvalue problem that 
can be solved iteratively. Even in the simplest cases 
considered by Hubbard & Brandenburg (2009), 7722 is 
proportional to 1/(1 — iwr), where r is the memory 
time. To understand this proportionality for large |w|, 
we note that for test fields (7) we find from (2) 

b' = iw“^ (iieqe„ -f 0(|a;|“^)) , 

where eq = k in the definition (48) of the operator 
Xleq. An illustrative example of the iterative proce¬ 
dure was given by Rheinhardt et al. (2014) for a more 
complicated case where a; = iA was complex. We can 
encounter a neutral dynamo such that ReA = 0 (this 
usually occurs for a specific value of fci) by increasing 
ki, i.e., decreasing the domain; see Figs. 1 and 2 of 
Rheinhardt et al. (2014) for a related problem. 

In the present case, because 7722 ( 271 ) is nonuniform. 


we have to allow for all possible wave numbers of the 
resulting mean field and compute the response for each 
wave number. This is just opposite to the usual mean- 
field dynamo problem and the MST approach where 
one computes the dynamo effects in the limit k —)■ 0. 
The relevant eigenvalue problem for our domain of size 
27r now becomes 


AA2{xi) 


iv + V22ixi,ki,iA))kl 

ki — 1 


X 



e*'=i(^i-«i)A(^i)dei. 


(61) 


In Fig. 16 we have already plotted 7722 ( 271 ,1,0), but we 
now need 7722 ( 271 , fci, iA) for all integer values of ki and 
a suitable value of w. 

Note that for our domain of size 27r the permissible 
wave numbers ki are integers. Furthermore, looking at 
the right panel of Fig. 15, we see that the eigenfunc¬ 
tion is odd about Xi = 0. This means that only odd 
values of ki contribute to the solution. In agreement 
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Figure 15 . The cca-averaged rms-normalised mean field 
b3{xi,X2)/{b3)^^‘^ from DNS for rj = 0.1 in a domain 
of size ( 27 r)^ (left panel) and the planar average &3(a:i) 
(obtained by averaging over X2 and X3, black line), 
and 53(3^2) (averaging over xi and a::3, red dashed line 
overplotted). Same colour-coding scheme as in Fig. 12 . 
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Figure 16 . r]22ixi) and 1733(3:1) for 77 = 0.1 and ki = 1 . 
7722(2:1) has a negative average, indicated by the dotted 
line. 


with our earlier experience, the amplitudes of the tur¬ 
bulent transport coefficients fall off quadratically with 
increasing wave number (see, e.g., Brandenburg et al. 
2008 b). We therefore expect that the compensated ex¬ 
pression 7722(2:1, fci, iA)fci should be independent of fci 
for large values. This is indeed the case, as can be seen 
from Fig. 17 , where we plot 7722(2:1, fci, 0 )A:i separately 
for odd and even values of ki. 

We should point out that these results are sensitive 
to the values of 77 and w, as will be demonstrated next. 
First, in Fig. 18 we plot 1722(2:1, A:i, iA) for w = 0 , 0 . 02 i 
and O.li, and for fci = 1 and 3 . While the shapes of 
the different curves remain similar, there is a signifi¬ 
cant reduction in the amplitude as ui increases. Thus, 
it is in general impossible to omit the memory effect. 
This agrees with earlier results for certain steady flows 
(Radler et al. 2011 ; Rheinhardt et al. 2014 ), although 
it is not a typical feature of turbulent flows (Hub¬ 
bard & Brandenburg 2009 ). Second, we give in Ta¬ 
ble 4 the volume-averaged values (7722) and (7733) for 
uj — 0 and different values of rj. It turns out that 
17 + (1722) = — 0.017 for 17 = 0 . 115 , tentatively sug¬ 
gesting that this case is weakly supercritical, while for 
77 = 0.120 we have 77 -I- (7722) = - 1 - 0 . 053 , which would 
be clearly subcritical. These values are close to those 
obtained from DNS, which show that the critical value 
of 77 for the onset of generation of magnetic field with 
the periodicities of the flow is around 0.1105 (see also 
Table 1 ). For more precise statements we would need 
to consider numerical solutions to ( 61 ). 

The fci-dependence of eddy diffusivity 77 -I- (7722) is 
shown in Fig. 19 and numerical values are given in Ta¬ 
ble 5 for 77 = 0.1 and a; = 0 . Here we compare the 
results from TFM obtained with the Pencil Code 
with those obtained by cancelling the exponential fac¬ 
tor and determining steady solutions to ( 59 ); the lat- 



Figure 17 . 7722(2:1, fci, 0 )fci for 77 = 0.1 shown colour- 
coded separately for odd (upper panel) and even (lower 
panel) values of fci. Same colour-coding scheme on 
both panels. 


ter were computed by the code by Fokkema ( 1995 ), 
employing the biconjugate gradients stabilised method 
BiCGstab(f) for £ = 6 (see Sleijpen & Fokkema 1993 ; 
Sleijpen & van der Vorst 1995 , 1996 ). We also show 
V + (7733) which is always positive. Note that rj + (7722) 
becomes zero at ki « 1 . 8 . Thus in our domain of size 
277 , where fci = 1 is the smallest wave number, the 
volume-averaged eddy diffusivity is clearly negative. 
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Figure 18. Compensated kernel rj 22 (xi, l,uj) (upper 
panel) and 7722 (a^i, 3,a;) (lower panel) for w = 0 (black 
solid line), 0.02i (red dotted line), and O.li (blue 
dashed line) using rj = 0 . 1 . 



Figure 19. Dependence of 77 +(?? 22 ) on ki, obtained with 
TFM by computing a steady solution to (59) (solid 
line), versus the results of the Pencil Code (filled 
symbols). For comparison, 77 + ( 7733 ) is shown (hollow 
symbols). Here, 77 = 0.1, a; = 0. 

6. Concluding remarks 

In mean-field electrodynamics, various analytical 
and numerical approaches are used to express the 
mean electromotive force, originally defined in terms 
of small-scale fluctuations of flow velocity and mag¬ 
netic field, as functions of the large-scale mean flow 
and magnetic field. Assessing the range of validity and 
clarifying conflicts in application of these approaches is 
crucial in view of many applications, e.g., in laboratory 


Table 4. ( 7722 ) and ( 7733 ) computed for w = 0 and vari¬ 
ous 77 . 


77 

{m 2 ) 

(^ 33 ) 

0.100 

-1.858 

0.230 

0.110 

-0.245 

0.239 

0.115 

-0.117 

0.242 

0.120 

-0.047 

0.244 


Table 5. Comparison between 77-|- (7722), obtained with 
TFM by computing a steady solution to (59) (first 
column, marked by an asterisk), and (7722) as well as 
( 7733 ) obtained with the Pencil Code for 77 = 0.1, 
a; = 0 and various ki. 


ki 

V + {m2)* 

{m2) 

{m3) 

2.5 

-fl.05656 

-fO.9559 

0.1843 

2.0 

-1-0.20545 

-f0.1055 

0.1843 

1.5 

-0.34576 

-0.4460 

0.1808 

1.0 

-1.75831 

-1.858 

0.2304 

0.5 

-2.11246 

-2.213 

0.3051 

0.25 

-2.02750 

-2.127 

0.3210 


experiments for dynamo generation or in astrophysics. 
For instance, a comparison of the traditional MFE ap¬ 
proach with those based on r-approximations of turbu¬ 
lence theory was carried out by Radler & Rheinhardt 
(2007). Here, we have compared two different meth¬ 
ods for estimation of the mean e.m.f.: MST, which 
explicitly considers steady or time-periodic laminar 
flows, and TFM, which is not affected by such a re¬ 
striction. For instance, Cabanes et al. (2014) recently 
observed in a rotating liquid sodium experiment “Der- 
viche Tourneur Sodium” a reduction of the effective 
magnetic diffusivity in some regions of the flow, prob¬ 
ably caused by the turbulence; it is thus important 
to assess, under which conditions the methods investi¬ 
gated here can be used to study this kind of problems. 

We have demonstrated that, in the two-scale setup, 
magnetic eddy diffusivities predicted by MST in three- 
dimensional small-scale steady flows are reproduced 
by TFM, provided volume averaging is used. It can 
be similarly shown that the same result holds true 
for time-periodic flows — in this case the averaging 
procedure must also involve time averaging over the 
temporal period. If other types of averaging (planar 
or over just one Cartesian variable) are applied, one 
can, in general, only expect a qualitative agreement 
between the results. One must also be aware of the fol¬ 
lowing caveat: Because of the asymptotic character of 
the MST results, achieving agreement with TFM eddy 
diffusivities requires small scale ratios e and thus high 
spatial resolution to be used when solving the TFM 
test problems (2). This can be seen as a drawback of 
TFM in comparison with MST. 

Results coinciding with those of MST can also be 
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obtained by a modified TFM algorithm that proceeds 
by setting a test field as initial condition and solv¬ 
ing the standard magnetic induction equation (60) 
without separating the field into mean and fluctuat¬ 
ing parts before the saturated regime for the magnetic 
field sets in, and then computing the e.m.f. due to the 
fluctuating flow and magnetic held as in the canonical 
TFM. Reliance of TFM on the integral (3) approxima¬ 
tion of the e.m.f., resulting in the ansatz (4) (or (57)) 
for the Fourier transforms, then proves crucial. This 
feature of TFM implies that, in the multiscale limit 
£ —)■ 0, the choice of the spatial variables for averaging 
becomes insignihcant, because the use of the Fourier 
transform over the remaining variables effectively con¬ 
verts all the averages into averages over all the three 
spatial variables (see (55) and (58)). This observation 
does not hold for the conventional version of MFE, 
where the e.m.f. is approximated by local differential 
operators. Nevertheless, the MST a-effect tensor in 
small-scale flows can be computed by TFM with the 
use of constant test fields and full spatial averaging. 

We have numerically confirmed the findings of Lan- 
otte et al. (1999) that the modified Taylor-Green flow 
possesses, in certain ranges of parameter values, nega¬ 
tive eddy diffusivities, and have shown that the same 
holds for the G.O. Roberts flow IV. This is in contrast 
with Devlen et al. (2013), who did not find negative 
eddy diffusivity for the former flow. Why did Lanotte 
et al. (1999) and Devlen et al. (2013) arrive to differ¬ 
ent conclusions for this flow? We have seen that the 
results of MST and TFM do not agree qualitatively un¬ 
less TFM applies volume averaging (R-IV depending 
on two spatial variables is a special case), but a num¬ 
ber of less important reasons make the picture even 
more complicated: (i) We have now obtained negative 
eddy diffusivities in mTG by TFM, but we have been 
forced to employ a planar averaging different from the 
one used by Devlen et al. (2013). (ii) Eddy diffu¬ 
sivities affecting the evolution of large-scale perturba¬ 
tions of distinct short-scale magnetic modes do not 
coincide. While Lanotte et al. (1999) considered eddy 
diffusivity for the neutral short-scale modes, Devlen 
et al. (2013) aimed at evaluating it for the dominant 
short-scale modes, which for rj = 0.02 is a distinct 
branch; hence, there are no reasons to expect their re¬ 
sults on eddy diffusivity to be interrelated. {Hi) Fur¬ 
thermore, at large scale separations (i.e., small e) the 
eddy diffusivity for the dominant branch for ry = 0.02 
is negative (see Fig. 6); in this case, TFM still can be 
used to evaluate the eddy diffusivity, but special pre¬ 
cautions must be taken in its implementation: TFM 
requires solving test problems (2), which are likely to 
inherit the instability of the unperturbed magnetic in¬ 
duction equation (60) and have exponentially growing 
modes. In the course of numerical integration, the 
growing modes will then set in due to the influence of 
round-off errors and progressively wipe out the contri¬ 


bution of the inhomogeneity in (2), which we intend to 
determine. A feasible strategy is to compute directly 
a time-independent solution to (2) regarded as a sys¬ 
tem of linear equations (after a suitable discretisation 
of the problem in space). 

Investigation of eddy diffusivity is supposed to yield 
the effective diffusivity which can be employed, e.g., 
to study nonlinear large-scale MHD regimes. Our re¬ 
sults demonstrate that this may be a non-realistic goal. 
Magnetic field in a nonlinear MHD regime can be de¬ 
composed into a linear combination of eigenmodes of 
the magnetic induction operator, where the coefficients 
are time-dependent. We have presented growth rates 
in branches of dominant large-scale magnetic eigen¬ 
modes generated by mTG in four symmetry subspaces. 
The parabolic shape of their plots as functions of 
the scale ratio £ near e = 0 or 1 confirms that the 
phenomenon of magnetic eddy diffusivity is observed 
for sufficiently large scale separations in all branches. 
However, Fig. 6 also shows that the curvature of the 
parabola varies significantly for different branches of 
modes. Hence, for a given molecular diffusivity, no 
universal eddy diffusivity tensor can be assigned to a 
given generating flow, because the action of eddy dif¬ 
fusivity significantly depends on the spectral compo¬ 
sition of the multiscale magnetic field itself, on which 
such an integral diffusivity acts. We are thus forced to 
conclude that a unified description of “average” mag¬ 
netic eddy diffusivity would only be possible in the case 
of a turbulent MHD regime with a well-defined statis¬ 
tics of the spectral composition of magnetic field — for 
instance, when a chaotic attractor of an MHD dynam¬ 
ical system is considered. (As a side remark, we note 
that the same holds true for the a-effect tensor: for 
large-scale permutations of different small-scale mag¬ 
netic modes generated by small-scale non-symmetric 
flows, generically different tensors are obtained, and 
hence the integral a-effect tensor depends on the spec¬ 
tral composition of the magnetic field.) From this 
prospective, the TFM approach seems advantageous, 
since it demands to separately quantify the influence 
of each Fourier harmonics in the mean magnetic field 
by nominating it as a test field and computing the fluc¬ 
tuating field that it induces and the respective output 
e.m.f., and afterwards to sum up such contributions of 
each individual harmonics into the integral output. 

Let us note some open questions, beginning with 
the most important one. 

1. TFM follows MFE in relying on approximations 
such as (3) (e.g., the precise kernels in (3) are not 
translation-invariant). It is essential, that neverthe¬ 
less in some limits (the most important of which is 
the limit of small magnetic Reynolds numbers) MFE 
yields results that are exact, i.e., TFM and MFE agree 
exactly and in all details with DNS (see the review by 
Brandenburg et al. 2010). Given that (i) we know 
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that MST is a precise corollary of the basic equations, 
and (ii) we have observed that MST can quantitatively 
disagree with TFM, we need to understand the math¬ 
ematical reasons for the aforementioned agreement in 
the respective limits. The availability of precise math¬ 
ematical demonstrations may help to determine the 
conditions, under which TFM can be used reliably. 

2. As we have mentioned in the Introduction, the 
generated large-scale structures are described by am¬ 
plitudes (depending exclusively on slow variables) of 
the small-scale neutral (magnetic or MHD stability) 
modes that constitute the leading term in the expan¬ 
sion of perturbation in the scale ratio. In two-scale sys¬ 
tems (such as the ones considered here), the temporal 
evolution of the amplitudes is governed by equations 
(mean-field or otherwise), where the a-effect opera¬ 
tor is never present together with the eddy diffusivity 
operator: since the orders of these differential opera¬ 
tors are different, they emerge at different orders of 
the scale ratio. (A joint action of molecular diffusivity 
and the a-effect is encountered in flows with an inter¬ 
nal spatial scale, see Chaps. 10 and 11 in Zheligovsky 
2011.) Do the two operators appear simultaneously in 
amplitude equations in a truly multiple-scale setup? 
In other words, can in such a setup the mean e.m.f. be 
a sum of the a- and vy-terms, as this was assumed in 
the early variants of MFE? 

3. Expressing entries of the magnetic eddy diffusiv¬ 
ity correction tensor in terms of solutions to auxiliary 
problems for the adjoint operator has proved useful not 
only for reducing the amount of computations, but also 
in analytical work, for establishing relations (29) be¬ 
tween the tensors for opposite flows, and for identifying 
zero entries of the tensor for translation-invariant flow 
(see Section 3). In these calculations we have relied on 
the similarity of the magnetic induction operator and 
the adjoint operator for the reverse flow. Do solutions 
to the auxiliary problem for the adjoint operator have 
a physical interpretation? 

4. Finally, the following technical question is of 
certain interest: Eq. (2) governing the evolution of the 
fluctuating magnetic field involves the operator 
where fp projects out the mean field and St is the usual 
operator of magnetic induction. Suppose that there 
is no small-scale dynamo action, i.e., all eigenvalues 
of St have non-positive real parts, and an averaging 
other than volume averaging is employed. Is it then 
possible, for some flows and some test fields, to have 
growing fluctuating solutions, i.e., can the operator 

have an eigenvalue with a positive real part? 
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